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Summary 

The three-dimensional, linearized Euler analysis, LINFLUX, is being developed to pro- 
vide a comprehensive and efficient unsteady aerodynamic scheme for predicting the aeroa- 
coustic and aeroelastic responses of axial-flow turbomachinery blading. LINFLUX couples a 
near-field, implicit, wave-split, finite-volume solution to far-field acoustic eigensolutions, to 
predict the aerodynamic responses of a blade row to prescribed structural and aerodynamic 
excitations. It is applied herein to predict the acoustic responses of a fan exit guide vane 
(FEGV) to rotor wake excitations. The intent is to demonstrate and assess the LINFLUX 
analysis via application to realistic wake/blade-row interactions. Numerical results are given 
for the unsteady pressure responses of the FEGV, including the modal pressure responses at 
inlet and exit. In addition, predictions for the modal and total acoustic power levels at the 
FEGV exit are compared with measurements. The present results indicate that the LIN- 
FLUX analysis should be useful in the aeroacoustic design process, and for understanding 
the three-dimensional flow physics relevant to blade-row noise generation and propagation. 
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1. Introduction 


The development of analyses to predict unsteady flows in turbomachines is motivated 
primarily by the need to predict the aeroelastic and aeroacoustic response characteristics 
of the blading. Accurate and efficient analyses are needed to determine the unsteady loads 
that act on the blades and the unsteady pressures that persist upstream and downstream 
of the blade row, for various sources of excitation. The latter include prescribed structural 
motions and/or aerodynamic disturbances at inlet and exit that carry energy towards the 
blade row. The computational resources required by nonlinear and viscous unsteady fluid- 
dynamic analyses continue to prohibit their use in detailed aeroelastic and aeroacoustic 
response studies. Therefore, approximate, e.g., linearized inviscid, analyses are needed to 
achieve efficient predictions of unsteady aerodynamic response phenomena. 

Much attention has been given in recent years to the development of three-dimensional 
(3D) linearized Euler analyses [HL93, HCL94, Sre96, MV97, MG98, MV98]. Such analyses 
account for the effects of nonuniform inlet and exit conditions, real blade geometry, steady 
blade loading, and transonic phenomena, on blade-row, unsteady aerodynamic response. In 
particular, the 3D LINFLUX analysis [MV97, MV98, VMC99] couples a near-field, implicit, 
finite-volume analysis for the unsteady perturbation of a general nonlinear mean flow in the 
near field, to analytic/numeric, far-field analyses for the unsteady perturbations of the mean 
flows, upstream and downstream, that depend only on radius. The near-field analysis is 
based upon the high-resolution, wave-split, finite-volume scheme [WJS88], employed in the 
time-dependent, nonlinear fluid-dynamic analysis, TURBO [Jan89, JHW92, CW93]. The 
far-field analyses, which are coupled to the near-field, finite-volume analysis at the near- 
field, computational inflow and outflow boundaries, allow incoming external aerodynamic 
excitations to be prescribed, and outgoing response disturbances to pass through the com- 
putational boundaries without reflection. 

A two-dimensional (2D) version of the LINFLUX analysis is reported in [MV95, VMK95] 
along with numerical results for subsonic and transonic unsteady flows through the 10th Stan- 
dard Cascade [FV93]. Comparisons of the 2D LINFLUX results with those of the potential- 
based linearization, LINFLO [VC84, UV91, HV91], and the nonlinear Euler/Navier-Stokes 
analysis, NPHASE [HSR91, SLH + 94], indicate that the former gives accurate response in- 
formation for the various types of unsteady excitation. 

The three-dimensional LINFLUX analysis is reported in [MV97, MV98, VMC99], and 
partially validated via result comparisons, with well-known 2D analyses, for benchmark 
unsteady flows. In particular, subsonic unsteady flows through a helical fan and a three- 
dimensional version of the 10th Standard Cascade (3D SC10), excited by prescribed blade 
motions, acoustic disturbances at inlet and exit, and vortical disturbances at inlet, have been 
studied. The helical fan and the 3D SC10 have high hub to tip ratios, and were subjected to 
highly two-dimensional unsteady excitations at low through moderate frequencies. Thus, the 
3D LINFLUX results could be compared meaningfully with those of the 2D classical analysis 
of Smith [Smi72] for the helical fan, and those of the 2D LINFLO analysis for the 3D SC10. 
The LINFLUX pressure responses, at midspan, are in very good agreement with those of 
the earlier 2D analyses. Moreover, these responses show reasonable radial trends, and the 
amplitudes and axial eigenvalues of the propagating and/or least-damped modal acoustic 
responses in the far field are accurately predicted. More recently [PVOO], LINFLUX results 
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for a standard computational aeroacoustics (CAA) configuration [Han99]; i.e., a flat-plate 
stator, with a hub-to-tip ratio of 0.5, excited by vortical excitations at inlet, were found to be 
in excellent agreement with the results of Namba’s [Nam87] 3D classical linearized analysis. 

The 3D LINFLUX analysis is being developed to serve as the unsteady aerodynamic 
component of an advanced design prediction system for fan noise. A fundamental source 
of such noise is the interaction between the wakes coming off the upstream rotor, or more 
generally, between the circumferential nonuniformities in the rotor exit flow, with the down- 
stream fan exit guide vane (FEGV). Current aeroacoustic design systems rely on unsteady 
aerodynamic analyses that are based on classical linearized theory, either 3D analyses or 2D 
analyses applied on radial strips. Classical analyses do not account for many of the impor- 
tant features of real flows. For example, the only stator configuration that can be addressed 
by a classical 3D analysis is one in which the stator consists of flat-plate blades that are 
aligned with a uniform mean or steady inlet flow. Because of such limitations, the unsteady 
aerodynamic element has been regarded as the weakest component in current design systems 
for predicting fan-noise generation and propagation. 

The LINFLUX analysis is being developed to overcome the limitations inherent to clas- 
sical linearized theory. In the present effort, we have applied this analysis to a realistic fan 
configuration to demonstrate and assess its capabilities for aeroacoustic design applications. 
We consider unsteady flows through the fan exit guide vane (FEGV) of a NASA/PW 22 inch 
advanced ducted propulsor (ADP) [Neu97, Pod97], operating at low (approach) and high 
(take-off) power settings. The unsteady flows through the FEGV are excited by circumferen- 
tial nonuniformities in the flow exiting the fan rotor, which are due primarily to the viscous 
wakes extending downstream from the rotor blades. The analytic wake excitation is based on 
measurements, taken downstream of the fan rotor, at NASA Glenn Research Center [Pod97]. 
A number of physical assumptions have also been introduced to provide an analytic wake 
excitation that is suitable for use with LINFLUX. LINFLUX is then applied to predict the 
unsteady aerodynamic responses of the FEGV to Fourier components, at multiples of the 
blade passing frequency (BPF), of the rotor wake excitations. 

It should be noted that the application of LINFLUX to predict the unsteady aerodynamic 
response of a blade row to a wake excitation involves the use four major computer codes: 
TIGER [SS91], to generate a 3D H-grid; AWAKEN [TE99, TopOO], to define analytic steady 
and unsteady inflow conditions based on measured or computed wake data; TURBO [Jan89, 
JHW92, CW93], to determine the steady background flow; and LINFLUX to determine the 
linearized unsteady flow. We will refer to this set of codes as the LINFLUX Prediction 
System, and provide instructions for running this system in [VerOO]. 

In this report, we describe the general unsteady aerodynamic problem for a blade row; the 
linearized unsteady aerodynamic equations, and the near- and far-field solution procedures 
that form the basis of the LINFLUX analysis; and the current wake-excitation model. We 
then proceed to give a detailed description of the application of LINFLUX to predict the 
rotor- wake/FEGV interactions in the NASA/PW 22 inch ADP. We will present numerical 
results for the steady and unsteady pressures and the unsteady velocity in the near field of 
the FEGV, and for the modal acoustic responses in the far field. Predicted results for the 
modal and total sound power levels at the FEGV exit will be compared with experimental 
measurements [Hei99]. 
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2. Unsteady Flow through a Blade Row 


We consider the general situation of time-dependent flow, of an inviscid perfect gas 
through a vibrating, rotor or stator, blade row (see Figure 1). The latter consists of N 
blades, which, in the case of a rotor, rotate about the duct axis at constant angular velocity 
and operates within a stationary annular duct with hub and duct, or tip, radii 
of r = r H and r = r D , respectively. In the absence of a vibratory motion, the blades are 
identical in shape, equally spaced around the rotor, and identical in orientation relative to 
an axisymmetric inlet flow. 

We will analyze this flow in a reference frame that is fixed to the blade row, in terms 
of cylindrical (r, 9 , £, t) and Cartesian (xi, X 2 , X 3 ,t) = (£, rsin0, —r cos 9, t) coordinates. 
Here £ and r measure distance in the axial flow direction and radially outward from the 
duct axis, respectively, 9 measures angular distance in the e ( xe r = e g direction, which, for 
a rotor, is opposite to the direction of blade rotation, and e is a unit vector. We will also 
have occasion to examine a rotor flow in the stator or absolute frame, and a stator flow in 
the rotor or relative frame. Thus, when needed for clarity, we will use the superscripts abs 
or rel to indicate that a physical quantity is measured relative to a stationary or a rotating 
reference frame, respectively. Thus, for example, 9 ahs = 9 Tel -f f }t. 

To describe flows in which the fluid domain deforms with time, it is useful to consider two 
sets of independent variables, say (x, t) and (x, t). The position vector x(x, t) =x + H(x, t) 
describes the instantaneous location, in the blade row frame of reference, of a moving field 
point, x refers to the reference or mean position of this point, and 7 Z(x, t ) is the displacement 
of the point from its reference position. The displacement field, 7£, is prescribed so that the 
solution domain deforms with the vibratory motions of the blades and is rigid far from the 
blade row. If the blades do not vibrate, is simply set equal to zero. 

In the present report, physical variables are, for the most part, given in dimensionless 
form. Lengths are scaled with respect to the reference length L Ref , time with respect to the 
ratio ^Ref/Fftef where Vj[ ef is the reference flow speed, velocity with respect to V^, density 
with respect to a reference density p Ref , pressure with respect to p^ ef (I/ R + ef ) 2 and specific 
internal energy with respect to (F Ref ) 2 . Here, the superscript * refers to a dimensional 
value of a flow variable. In LINFLUX, the reference length is taken to be the blade chord 
at midspan; the reference fluid density and flow speed, to be the inlet, temporally- and 
circumferentially-averaged density and relative (to the blade row) flow speed, at midspan, 
respectively. 

For aeroelastic and aeroacoustic applications, we are usually interested in a restricted class 
of unsteady flows; those in which the unsteady fluctuations can be regarded as perturbations 
of a background flow that is steady in a reference frame that rotates with the blade row. 
Moreover, the steady background flows far upstream (£ < £_) and far downstream (£ > £ + ) 
from the blade row can be assumed to consist of at most small steady perturbations from 
steady flows that depend only on radius. The time-dependent or unsteady fluctuations 
in these flows arise from temporally and circumferentially periodic unsteady excitations 
of small-amplitude, i.e., prescribed vibratory blade motions and prescribed aerodynamic 
disturbances at inlet and exit that travel towards the blade row. 
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For example, if the blades vibrate at reduced frequency, u, as seen by an observer in the 
blade-row frame, and at constant interblade phase angle, a, we can write 

0 + 27tti/N, £, t) — T n i?e{R B (f, 9, £) exp[i(utf + na)]}, x on B . (2.1) 

Here, 7l Bn is the displacement of a point on the nth moving blade surface from its mean 
position in the rotating frame; T n is a rotation matrix, which relates a vector in the reference 
(n = 0) passage to its counterpart in the nth passage; n = 0, 1,2 , ...,N - 1 is a blade 
index; Re{ } denotes the real part of { }; R B is the complex amplitude of the reference 
(n = 0) blade displacement; and B refers to the mean position of the reference blade. 
The interblade phase angle, cr, can be determined from the nodal diameter pattern of the 
unsteady disturbance. Thus, for the vibratory blade motion, described in equation (2.1), 
<7 = 27tNd/N, where |A^d|, is the number of nodal diameters, i.e., the the number of times 
the vibratory disturbance repeats around the wheel. The sign of N D depends on the direction 
of rotation of the disturbance pattern. If the vibratory disturbance travels in the negative 
^-direction, Np > 0. Thus, for a rotor, if the vibratory disturbance travels in the direction 
of blade rotation, No > 0. 

The unsteady disturbances in the far upstream and far downstream regions are, in part, 
prescribed as a fluid dynamic excitation and, in part, depend upon the interaction between 
the fluid and the blading. Typically, an unsteady aerodynamic excitation is represented 
as a linear combination of fundamental disturbances. The fundamental disturbances are 
harmonic in time, at temporal frequency u, and circumferentially, at angular wave number 
fh = iV D + mN , m = 0, ±1,±2, . .. . For example, if the underlying absolute mean or 
steady flow is uniform, the pressure associated with a fundamental acoustic excitation is of 
the form 

Pa( t > e > £> *) = Rei^p^ir) expfc^ + i(md + wf)]}, £ > . (2.2) 

Here, p^{x.,t) describes an incident pressure disturbance, i.e., a pressure disturbance that 
travels towards the blade row from far upstream (£ < £_) or far downstream (£ > £ + ). The 
quantities w, m = iV D + mN or a m — 2 tt(N d /N + to), and the disturbance amplitude, 
a ± , are prescribed; the radial mode shape, p R,± (r), and the axial exponential coefficient 
X ± = P ± + i K f > where /3 is the axial attenuation coefficient and is the axial wave number, 
are determined from the equations that govern the unsteady fluid motion in the far field. 
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3. Unsteady Aerodynamic Model 


Nonlinear field equations that govern the unsteady flows, described in §2, are derived from 
the mass, momentum and energy conservation laws for inviscid flow and the thermodynamic 
relations for a perfect gas. These equations, when referenced to a coordinate frame that 
rotates with the blade row, at constant angular velocity S7, can be written in integral form 

^ . f r r 

- / U dV+ [F*. - \3il Xj )n Xj dA = / SdV . (3.1) 

dt Jy J A Jy 

Here, the symbol ~ indicates a time-dependent flow quantity, and V(t) is an arbitrary 
moving control volume, which is bounded by the control surface -4(x, t) = 0. Also, 7?.(x, t) = 
{7l Xl , 7l X2 , Ti X3 ) is the prescribed velocity of the moving field point at x = x + 7?.(x, t), n is 
a unit vector normal to the surface, A, and pointing outward from V, and a summation over 
repeated indices is implied. The source term on the right-hand-side of (3.1) accounts for the 
rotation of the reference coordinate frame. 

The state, U, flux, F x ,j — 1, 2, 3, and source term, S, vectors in equation (3.1) are given 
by 


' _~p 


v j+ 1 1 


0 

pV Xl 


U j+l U 2 fUi + PSij 


0 

Wx2 

. F*,(U) = 

U j+l Uz/Ui + PS 2j 

, S(U,x) = 

Q 2 U x x 2 4- 2f IU 4 

pv x3 


Uj+lUi/Ui + P5 3j 


n 2 i \x z - 2 nu 3 

. P&T . 


. U j+1 (U 5 + P)/U x . 


. n 2 (u 3 x 2 + u 4 x 3 ) . 


(3.2) 

Here, V, Er = E + V 2 /2 and P = ( 7 - 1)1) E = (7 - l )[£/ 5 - + U'i)/2} 

are the time-dependent fluid density, velocity, specific total internal energy, and pressure, 
respectively. The velocity and total internal energy are measured relative to the blade-row 
frame of reference. 

We also require a differential form of the field equation (3.1), which applies at fixed 
locations ('R. = 0 ) in the blade-row frame, to describe the unsteady flows far upstream 
(£ < f_) and far downstream (£ > £+) of the blade row. Expressed in terms of cylindrical 
coordinates, this equation has the form 


dU cyl 

dt 


+ r 


-A* +r -.*& + «£«S. 


dr 


dd 


(3.3) 


The state and source-term vectors in (3.3) are given by 


U cyl = < 


~p ] 


' 0 


r 0 

pV r 


(^ yl ) 2 /C^ yI + P 


2 Uf + Ul yl Qr 

pVg 

> , S = r 1 < 

-ufuf/uf 

> -(- < 

-2Uf > 

~pVj 


0 


0 

, P&T j 


0 J 


UfClr- J 


(3.4) 
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The flux vectors F r (U cyl ), FV(U cyl ) and F^(U cyl ) and the pressure P(U cyl ) have functional 
forms similar to those indicated previously for F Xj .(U), j = 1,2,3, and P(U). 

The foregoing field equations must be supplemented by conditions on the unsteady flow 
at the blade surfaces, the duct walls, and the inflow (£ = £_) and outflow (£ = £ + ) boundaries 
of a near-field computational domain. Flow tangency conditions, i.e., 

(V — 7V) • n = 0 for x € B n , r — r B and r = r D , (3.5) 

apply at the moving blade surfaces, B n , and at the stationary {it = 0) duct walls. In 
addition, temporally- and circumferentially-averaged values of the total pressure, the total 
temperature and the inlet flow angle are specified as functions of radius at the computa- 
tional inflow boundary, and the circumferentially- and temporally-averaged static pressure is 
specified at the outflow boundary, consistent with radial equilibrium. Unsteady departures 
from these averaged flow conditions at inlet and exit that carry energy towards the blade 
row must be specified; those that carry energy away from the blade row must be determined 
as part of the solution for the unsteady flow. 

3.1 Linearized Unsteady Aerodynamic Equations 

Since the unsteady excitations are assumed to be of small amplitude (e.g., \R, Bn \ ~ 
@( e ) « 1) and to occur at a single temporal frequency, u, the unsteady component of the 
inviscid flow can be approximated as a first-order (in e), or linear, harmonic perturbation 
of an underlying zeroth-order, nonlinear, background flow. The background flow is steady 
in the blade-row frame. This approximation offers several computational advantages. First, 
since the first-order unsteady flow properties have harmonic time-dependence, physical time 
dependence, can be removed from the resulting unsteady boundary-value problem. Second, 
as a consequence of our assumptions regarding rotor geometry, the inlet and exit mean- 
flow conditions, and the circumferential behaviors of the unsteady excitations, numerical 
resolutions of the steady and linearized unsteady flows can be limited to a single, extended, 
blade-passage region; i.e., a region of angular pitch A 0 - 2n/N. Finally, the steady flows far 
from the blade row, depend only on radius. Therefore, approximate representations of the 
unsteady perturbations in the far field can be constructed, and matched to a computational 
near-field solution, to limit the axial extent of the near-field, computational domain. 

To determine the zeroth- and first-order aerodynamic equations, we expand the unsteady 
state vector, U, into an asymptotic series of the form [HC93] 

U[x(x,t), f] = U(x) + u[x(x, *),<] + — = U(x) + Re{ u(x) exploit)} + . . . . (3.6) 

The column vectors U(x) and u[x(x, t),t] contain the conservation variables for the steady 
background flow at the mean position, x, of a moving field point and the first-order unsteady 
flow at the instantaneous position, x = x + Tt{x, t) = x + Re{R(x) exp(fo;t)}, of this point, 
respectively. The components of the vector u are the complex amplitudes of the first-order 
unsteady conservation variables, i.e., u T = [p , pv Xl +pV Xl , pv X2 +pV X2 , pv X3 + pV X3 , pe T + 
pE T ] where p, V and E T and p, v, and er are the steady and the complex amplitudes of 
the first-order, unsteady, primitive, flow variables. The unsteady flux F Xj and source term, 
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S, vectors are expanded into Taylor series about the mean flow state, U, and the reference 
spatial location, x, i.e., 

dF aq 

F*,(U) = F,,(U) + -gJiu+... and S(U,x) = S(U,x) + ^u+(W- V s )S+ . . . . (3.7) 

The integral field equation (3.1) is expressed in terms of the moving control volume, V, 
and the moving control surface, A. The corresponding differential field equation would be 
expressed in terms of the moving spatial coordinate x. However, because of the dependent 
variable expansion (3.6), it is more appropriate to express the steady and linearized unsteady 
equations in terms of the mean or steady-state quantities, V, A, and x. To within first-order 
in e, the required spatial transformation relations are 


dV — (1 + Vx • 7 Z)dV + . . . , ndA = noL4 -I- A(noL4) exp(zud) + . . . 
and 

^ /dxj = d fdxj - (dn Xm /dx,)d /dx m + . . . , 


(3.8) 


where (^7x • Tty dTZ m /dx m , n is the unit outward normal vector to the control surface 
A, and A(ncL4) is the complex amplitude of the first-harmonic component of ndA — ndA. 
Finally, to within first order in e, the temporal derivatives appearing in the nonlinear, time- 
dependent, integral and differential field equations transform according to 


d 

dt ,/y 


/,< 


fd () 

+ ( ) Vx • 71 ) dV + . . 

and ~ 

d 

\ dt 

s / 

dt 

x ~ dt 


+ Ttxr, 


d 

dx„ 


+ 


(3.9) 

The equations that govern the zeroth-order steady and the first-order unsteady flows are 
obtained by substituting the foregoing series expansions into the time-dependent, nonlinear, 
governing equations; equating terms of like power in e; and neglecting terms of second and 
higher order in e. This procedure leads to nonlinear and linear variable-coefficient equations, 
respectively, for the zeroth- and first-order flows. 

The nonlinear, integral, field equation that governs the steady background flow is 



F Xj T^xj dA 



In addition, the flow tangency conditions, 


(3.10) 


V • n = 0 , for x 6 B n , r = r H and r = r D (3.11) 

apply at the mean blade surfaces and the duct walls, and circumferential periodicity condi- 
tions; e.g., 

p(f,9 + 2nn/N,O = p{r,0,O and V(f, 6 + 2im/N, l) = T n V(f, 6, £) (3.12) 

apply upstream and downstream of the blade row. Finally, circumferentially averaged values 
of the appropriate steady flow variables are specified as functions of radius at the inflow and 
outflow boundaries and circumferential harmonics of these variables are allowed to evolve to 
values that are consistent with a blade row operating within a long annular duct. 
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The linear equation that governs the first-harmonic unsteady perturbation is 

iw i udV + L ( w Lu ) n “ iA - L £? ad? - - L " (Hr ) vdV 

(3l3) 

+ / n Xj dA - F,,A(n<M)] + / JL(r S)dV , 

J A Jv vXj 

where the variable coefficients depend on the steady background flow, and the source terms 
depend upon this flow and the prescribed displacement field, R. 

Linearized flow tangency conditions, i.e., 


v • n = iu R • n + V • V(R • n), for x € B n , r = r H and r = r D , (3.14) 

apply at the mean blade surfaces and at the duct walls, and phase-shifted, circumferential, 
periodicity conditions, e.g., 


p{r, 9 + 2-irn/N , f) = p(r, 0 , 0 and v(f, 0 + 2nn/N, = T n v(f, 9 , f) exp(mcr) , (3.15) 

apply upstream and downstream of the blade row. The far-field conditions imposed in 
the unsteady problem must allow for the prescription of external unsteady aerodynamic 
excitations and permit unsteady disturbances coming from within the solution domain to 
pass through the computational inlet (at £ = £_) and exit (at £ = £ + ) boundaries without 
reflection. 

Differential field equations that describe the steady and the first-order unsteady flows at 
fixed points, x = x, in the far upstream (£ < £_) or far downstream (£ > £ + ) regions, can be 
determined by applying the series expansions (3.6) and (3.7), with R = 0, to the nonlinear, 
time-dependent equation (3.3). We find that 


and 


, drF r 

r _ 1 -r— + r- 

or 


_id(r Au) 
luju + r — hr' 


dr 


dFg dFs 

00 ae 

(3-16) 

0Bu 0Cu _ 


d$ + 0£ ° U ~ ° ’ 

(3.17) 


where A — £?F r /dU cyl , B = dF$ /9U cyl and C = /d\J cyl are flux Jacobian matrices 
and D = (9S /<9U cyl is the source-term Jacobian. 


3.2 Solution Strategy 

We require sequential solutions to the foregoing nonlinear steady and linear unsteady, 
boundary value problems, over a single, extended, blade-passage, solution domain, to predict 
the unsteady aerodynamic response of a blade row to a prescribed unsteady excitation. In the 
present study, we apply the nonlinear analysis, TURBO [JHW92], to determine the steady 
background flow, and then seek solutions to the linear unsteady problem by matching a wave- 
split, finite- volume analysis for the unsteady perturbation in the near field, i.e., in the region 
£-<£<£+, to approximate analytic/numeric descriptions of the unsteady perturbations 
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of the radially-dependent , steady background flows that exist in the regions far upstream 
(£ < £-) and far downstream (f > £+) of the blade row. 

For the near-field, finite- volume analysis, it is advantageous to regard the state vector u 
as pseudo time dependent, i.e., to set u = u(x, r), where r is the pseudo time variable, and 
add the term d(f 9 u dV)/dr to the left-hand side of equation (3.13). This allows the use of a 
conventional time-marching algorithm to converge the solutions for the complex amplitudes 
of the first-order, unsteady, flow quantities, to constant values. 

The complex amplitude of the displacement field, R(x), must be prescribed over the entire 
solution domain. This field is defined so that the solution domain deforms with the blade 
motion (i.e., R(x) = R Bn (x) for x € B n ), slides along the hub and duct walls (R(x) • e r = 0 
for f = r H , r D ), and remains rigid far from the blade row (R(x) = 0 for £ > In addition, 
R(x) is prescribed along one blade-to-blade periodic boundary, and set at the other boundary, 
so as to satisfy phase-lagged periodicity, cf. (3.15). In the near field, R(x) is determined, 
first, along the hub and duct walls, and then, in the interior of the computational domain 
as solutions of Laplace’s equation, V?R = 0, subject to the appropriate Dirichlet boundary 
conditions. For unsteady flows in which the blades are stationary in the blade-row frame, R 
is simply set equal to zero. 
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4. Near-Field Numerical Model 


The numerical procedures used in LINFLUX to resolve the linearized unsteady flow 
in the near field, are based on those employed in the nonlinear analysis, TURBO [Jan89, 
JHW92, CW93]. Both analyses use a sheared H-mesh [SS91], which defines a curvilinear 
coordinate system, such that, there is a one-to-one correspondence between the points, x, in 
the physical domain and the points, a, in a rectangular computational domain, where the 
grid is uniform and orthogonal. The q x , q 2 and a 3 computational coordinates, or the /, J, K 
computational mesh indices, correspond, generally, to the axial, radial and circumferential 
directions, respectively. Cell faces are surfaces of constant computational coordinate, with 
each cell bounded by the six surfaces: a \ = / ^ 1/2, « 2 = J T 1/2 and a 3 = K =F 1/2. 

4.1 Finite Volume Equations 

If we let the symbol * refer to a quantity expressed in terms of cell geometry, then a finite- 
volume spatial discretization of the pseudo time-dependent form of the linearized unsteady 
equation (3.13) can be written as 

Ag 

~ = —iujix - Sjfj + — u - zw(Ai?)U - 5j{f -I- (Atf)S + $(R • V*)S = -r , (4.1) 

where u = t?u, f, = (dFj/dU)\x = A jk (d F k /dXJ)u = A jk i k , and if = -tfjU + a jk F k . Here, 
the quantities u, U, and S represent average values of the physical state and source term 
vectors over a mean cell volume; d is the mean cell volume; Aj k is the mean area of a constant 
otj cell face projected in the x k direction; the vectors Fj and fj are the steady and unsteady 
fluxes, respectively, across a constant otj cell face; if is the unsteady flux associated with the 
displacement field, R; and r is the residual of the first-harmonic unsteady equation. The 
steady quantities d, Aj k , U, Fj(U) and S(U, x) and the unsteady displacement field, R, 
are regarded as known for the linearized unsteady analysis. 

The operator 5j in equation (4.1) denotes the difference in the /-direction across adjacent 
cell interfaces, and the repeated j index implies summation over all computational coordinate 
directions. Thus, the terms Sjij and Sjif are the net unsteady fluxes through a cell due to 
the unsteady fluid motion, and the grid motion, respectively. The linearized perturbation 
equation contains source terms that are associated with changes in cell volume, cell face 
area, and cell radial location. These terms depend on known steady flow properties and 
the prescribed displacement field, R(x). In evaluating the volume terms At? = Sj(A jk R Xk ) 
and t?j = iu)Aj k R Xk , the R Xk are taken to be the average displacements over a cell face. 
The perturbations in the cell face areas, a,j k , and in cell location are computed using the 
displacements of the cell vertices. 

The unsteady field equation must be solved subject to conditions at the boundaries of the 
near-field computational domain. The flow tangency conditions at the blade surfaces and 
duct walls, cf. (3.14), are implemented using phantom cells inside a solid surface. The den- 
sity and pressure in a phantom cell are defined by a first-order accurate reflection condition, 
and the phantom cell velocity is defined such that the velocity at the solid surface, which 
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is the average of the velocities in the phantom and the interior cells, satisfies the flow tan- 
gency condition, in a manner consistent with the finite volume discretization. Phase-shifted 
periodicity, cf. (3.12), is imposed along a passage boundary by taking advantage of grid pe- 
riodicity in the blade-to-blade direction. Thus, we set u| 7 JK = T±iUi y j t jc±(NK-i) exp(±«a) 
where NK denotes the number of blade-to-blade grid points, and the T„ matrix rotates 
the x 2 and x 3 components of the momentum vector through n blade passages. Finally, 
analytic/numeric far-field solutions are matched to the numerical near-field solution at the 
computational inflow and outflow boundaries, i.e., at £ = 

After replacing the pseudo time derivative in (4.1) by a first-order accurate, two-point, 
backward, difference expression, and expanding the residual, r n+1 , about the nth time level, 
we can write the discretized unsteady field equation as 

(ei + Jj) Au" = (e + iu$ - .»H) Au" + S, (fb A u"j = — r» . (4.2) 

In equation (4.2), the superscript n refers to the nth pseudo time level, © = d/A t, Au n = 
u n+1 — u", f is defined as in (4.1), and dr/du is a constant, since the unsteady residual is a 
linear function of the state vector u. 

4.2 Evaluation of Flux Terms and Pseudo-Time Integration 

To simplify the descriptions of the spatial discretizations, used to approximate the flux 
terms in equation (4.2), and the pseudo-time integration, used to solve this equation, we will 
consider a “one-dimensional” flow in which Fj = F and f) = f are the steady and unsteady 
flux vectors in the aj = a coordinate direction. The subscript J refers to the cell volume 
bounded by the cell surfaces ata = J-f-l/2 and a = J — 1/2. 

Flux difference splitting is applied to evaluate the flux terms in equation (4.2). The flux, 
fj+i/ 2 > at the J + 1/2 cell interface is evaluated in terms of the flux in the cell to the left (J) 
of the interface and the flux due to waves approaching the interface from the right. Thus, 

(u/ +1 - u j) , (4.3) 

uj?! 1/jA+i/i 

where i(nj,A J+ \j 2 ) is a flux based on the state vector in the Jth cell and the mean area 
of the J -I- 1/2 cell interface, and the (— ) superscript on dFjd\3 indicates left traveling 
disturbances, which correspond to the negative eigenvalues of the flux Jacobian matrix. The 
flux Jacobian matrix d¥ /<9U|uRoe^ ^ is evaluated in terms of the Roe-averaged [Roe81] 

intermediate state vector, / 2 , and the area, A J+1 / 2 . 

The discrete approximation (4.3) is first-order accurate, since the interfacial fluxes are 
based only upon information from adjacent cells. Second order spatial accuracy is achieved 
by adding corrective fluxes, which bring in information from additional neighboring cells. In 
LINFLUX, the corrective perturbation flux at the J + 1/2 interface is comprised of right 
traveling waves at the upstream interface ( J — 1/2) of the Jth cell and left traveling waves 


V+l/2 


dF 

dU 


u 


J +1/2 


f(uj, Aj +1 / 2 ) + — 
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at the downstream interface (J + 3/2) of the ( J + l)th cell. These waves are approximated 
using the Roe-averaged flux Jacobian matrix at the J + 1/2 interface. 

Once the interfacial fluxes have been computed, they are spatially differenced; i.e., 

5i ss ij+\/2 - fj— 1 / 2 > to compute the net flux terms that appear on the right- and left-hand 

sides of equation (4.2). A second-order discrete approximation is used to evaluate the net 
unsteady flux term that appears on the right-hand side of (4.2). The left-hand side flux 
term, <5[(dF/dU)Au] is evaluated using the first-order accurate approximation in (4.3). 

The linearized unsteady equation (4.2) is discretized, as outlined above, leading to a 
system of linear algebraic equations. This system is solved, at each pseudo-time step, using 
a symmetric Gauss-Seidel (SGS) iteration procedure in which the left-hand side matrix is 
decomposed into diagonal and off-diagonal submatrices. Thus, we set 

(© jI +|^ j All} = DjAu" - M+_! Au}_ x + M7 +1 Au} + 1 = -rj , (4.4) 

where the D submatrix contains the diagonal elements of the original matrix, and the M + and 
M - submatrices contain the off-diagonal elements in the negative and positive computational 
coordinate directions, respectively. 

Solutions to equation (4.4), which are coupled iteratively to analytic/numeric far-field 
solutions, are marched forward in pseudo-time until convergence is achieved, i.e., ||r n || -+ 0. 
The current LINFLUX implementation uses explicit surface boundary conditions. These are 
incorporated into the SGS iteration procedure, so that they are imposed in a semi-implicit 
manner. This treatment has been found to yield better convergence properties than a purely 
explicit implementation. 
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5. Unsteady Perturbations in the Far-Field 


The unsteady solution must be extended into the regions far upstream (£<£_) and far 
downstream (£ > £ + ) from the blade row. To meet this requirement, we will match approx- 
imate far-field solutions, determined from a reduced set of linearized unsteady equations, 
to the near-field finite-volume solution. The far-field solutions enable the prescription of 
unsteady aerodynamic disturbances that enter the computational domain (excitations), and 
the determination of the outgoing disturbances, including the outward propagating acoustic 
response disturbances. 

We have assumed that, far from the blade row (£ > ^ T ), the inner, and outer, r D , 
duct radii are constant, and the circumferentially-averaged, steady flow depends only on 
radius. Thus, it follows from the governing equations, that the steady velocity has axial 
and circumferential components, but no radial component, and dP/dr = p{Vg hs ) 2 /r. Under 
these conditions, the linearized unsteady field equation (3.17) reduces to 


(5.1) 


Equation (5.1) describes the first-order unsteady perturbations far from the blade row. The 
subscript 2 on the Jacobian matrices in (5.1) indicates that they are evaluated at U £ yl = 
pV r = 0; e.g., A 2 = dF r /dU c %cy. =0 . 

If the absolute mean flow is uniform in the far field, analytic solutions to equation (5.1) can 
be determined [TS62, VTM82, VMK95]. Such solutions indicate that an arbitrary first-order 
unsteady aerodynamic perturbation can be represented as the sum of independent entropic, 
vortical and irrotational acoustic disturbances. The entropic and vortical disturbances are 
convected by the mean flow, and, if the axial mean flow is subsonic, modal acoustic dis- 
turbances propagate upstream and downstream, or attenuate with increasing axial distance 
from the blade row. 

For nonuniform mean flows, the situation is much more complicated [Kou95, GA97, 
TA98]. In particular, for rotational, but isentropic, mean flows, the unsteady entropy is an 
independent convected disturbance, but the unsteady vorticity and pressure are coupled. As 
a result, nearly-convected, vorticity-dominated, modal disturbances, that contain pressure, 
and acoustic or pressure-dominated modal disturbances, that contain vorticity, occur. For 
nonisentropic and rotational mean flows no purely convected disturbance field will exist, 
because the entropic and vortical disturbances are coupled to the acoustic disturbances 
[Ker77]. 

It is difficult to determine the axial and radial behaviors of nearly-convected modal 
disturbances and very difficult to represent an arbitrary far-field disturbance in terms of 
series of acoustic and nearly convected modal disturbances. Therefore, to enable practical, 
three-dimensional, unsteady-flow calculations, some liberties must be taken in describing the 
unsteady perturbations in the far field. 

In the present effort, these perturbations will be modeled as the sum of a finite series of 
acoustically-dominated modal disturbances and a “convected” or “rotational” disturbance 
field. Thus, we set u = 4- Uc, where the state vector describes the modal acoustic 
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disturbances, and describes the vortically and/or entropically dominated unsteady per- 
turbation. The state vector u c is prescribed at inlet (£ = £_) as an approximate solution to 
(5.1), and extrapolated or convected through the exit boundary (£ = f + ) as a solution of the 
convection equation Duc/Dt — 0. The axial and radial behaviors of the modal acoustic dis- 
turbances are determined by far-field eigenanalyses. The amplitudes of the incoming acoustic 
disturbances are prescribed; those of the outgoing acoustic disturbances are determined by 
coupling the near- and far-field solutions for the unsteady state vector u. 

5.1 Acoustic Eigensolutions 

The acoustic disturbance is determined by substituting an assumed form for a modal 
disturbance, i.e., u = u* n (r) exp(x mn £ + ifh8), into the field equation (5.1), leading to a 
system of homogeneous ordinary differential equations of the form 

d 

iul u« n + t~ 1 — (rA 2 ul) + iffir- 1 B 2 u£ n + Xmn C 2 - D 2 u£ n = 0 . (5.2) 

After discretizing this system, and applying the duct wall boundary conditions v r = 0 at 
t = rjj and r = r^, we obtain the matrix equation 


(P “ XmnC 2 ) Ul = 0, (5.3) 

where P = — iul — L(r, A 2 ) — imr~ l B 2 +D 2 and L(r, A 2 ) is a finite difference approximation 
to r~ l d(rA 2 u^ n )/dr. The column vector u^ n in equation (5.3) contains an entry for each 
of the five conservation variables at each radial discretization point. 

Equation (5.3) can be solved, using a standard linear algebra routine, to determine 
the axial eigenvalues, Xmn, and the associated right eigenvectors, u^ n (r), of the modal, 
acoustic disturbances in the far-field. The left eigenvectors, u^ n , are determined by solv- 
ing the equation (P - Xmn C 2 ) H u^ n = 0, where the superscript H denotes the conju- 
gate transpose. An orthonormal set of left eigenvectors can then be obtained by setting 

(' V mn) H = (Umn)*C 2 /[( u£j"C 2 u£„]. 

The numerical solutions to (5.3) may contain spurious modes, which must be eliminated 
to yield a valid solution set. In addition, both acoustic and nearly convected disturbances 
can emerge as solutions to (5.3). The group velocity 


V - duJ - ( V mn 1 Cgugn) ,, 

9 ' mn d Xmn (y^idP/d^nU ’ (5 - 4) 

i.e., the axial velocity at which an unattenuated mnth modal disturbance carries energy, 
can be used to distinguish between such solutions. For example, nearly-convected distur- 
bances travel downstream, without attenuation, at axial speeds slightly less than and slightly 
greater than the mean flow speed. If the axial mean flow is subsonic, acoustic disturbances 
travel both upstream and downstream. Since the modal nearly-convected disturbances are 
presumed to be included within the disturbance field represented by the state vector u c , 
they must be eliminated to determine a valid acoustic solution. 

After the spurious and nearly-convected modes have been eliminated, the remaining 
eigen-disturbances can be combined to construct the acoustic state vector, u^. The latter 
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can be written in the form 

M M 

u A (r, e, 0 = Y ex P(^ 0 ) Y u ( r ) ex P + a rn^ (r) exp (xiOl 1 

(5.5) 

to emphasize the distinction between upstream, denoted by the superscript - , and down- 
stream, denoted by the superscript 4- , traveling waves. Here, M and Mu are positive 
integers; the index fx indicates the number of radial nodes or zero crossings in the fxth radial 
mode; the are the modal complex amplitudes; the u*’* (r) are right eigenvectors; and 
the x*. = ^ are axial eigenvaiues. 

The radial modes in (5.5) are ordered by their number of zero crossings. To determine 
if any of the lower order modes are missing, or if spurious modes have been kept, one can 
count the zero crossings in the kept modes. If there is one upstream and one downstream 
traveling mode for a each number of zero crossings, then it is likely that only the genuine 
modes have been retained. 

The complex amplitudes, a of the outgoing modal disturbances are determined by 
taking inner products involving the and the near field state vector u; i.e., 



Thus, it is being assumed that the inner products (v^, u c ) are negligible; i.e., that the left 
eigenvectors of the acoustic disturbances are nearly orthogonal to the convected disturbance 
field. 

The far field solutions must be applied in conjunction with a numerical near-field solution 
to determine the linearized unsteady flow. In the near field, the linearized unsteady governing 
equations are solved using the pseudo-time marching technique described in §4. After each 
iterative update of the near-field solution, the amplitudes of the outgoing acoustic modes, 
i.e., a~ at £ = and at £ = £+, are determined by taking inner products, cf. (5.6), 
using the near-field state vector, u. At the downstream far-field boundary, the acoustic 
disturbance, u^, is subtracted from the total unsteady disturbance, u and the remainder, 
u c — u — ux, is regarded as a convected disturbance; i.e., as a solution of the equation 
Duc/Dt = 0, to provide an update for Uc(r, 6, £+) at the outflow boundary. The new far- 
field solutions, which are based on the finite numbers of circumferential and radial acoustic 
modes that can be represented accurately on the near-field computational grid, are then used 
to supply the inflow and outflow information for the next near-field update. 

5.2 Acoustic Response Parameters 

It is desirable to introduce real parameters that provide useful descriptions of the un- 
steady aerodynamic response. For example, the work per cycle, which describes the work 
done by an airstream, on a blade, over one cycle of motion, is an important parameter in 
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blade flutter investigations. For aeroacoustic applications, the sound pressure and sound 
power levels provide convenient measures of the acoustic response far from a blade row. 
Expressions for evaluating these quantities are derived below. 

Consider the first-order unsteady flow variables, 5(x, t) = Re{a(x) exp(ia>t)} and 6(x, t ) = 
Re{6(x) exp(iud)}. The time-mean value, (56), of the product 56, is given by 

r2i r 

{ab) = (27T)- 1 / abd{ut) = {ab* + a*b )/ 4 = Re{o672} , (5.7) 

JO 

where the superscript * indicates the complex conjugate. The area-averaged, time-mean 
value, ((ab)}, of ab is given by 

((ab)) = j (ab)dA D = ^Ap l Rel^Jr j ab*rd6dr^ . (5.8) 

If the flow variables d and 6 are associated with the far-field acoustic response of a blade 
row, then, it follows from (5.5) that the product ab* can be expressed in the form 


OO CO OO OO 

ab* = 5Z ^ a ^ a n^m^nuY x exp [i(m - n)fl + (Xm|. + XwK] » ( 5 - 9 ) 

T 7 X— — OO (i= 0 n=-00 V— 0 

where a mfl and a nu are modal complex amplitudes and ^^(r) and ^\ v (r) are modal ra- 
dial shape functions. After substituting (5.9) into (5.8), integrating with respect to 6, and 
restricting consideration to propagating acoustic response waves, we find that 


{{ab}) = nApRe 



Q'miiQ'rnv 


exp [i(K mii - «„,„)£] x /: rmMlSrdv 1 . (5.10) 


For acoustic perturbations of uniform absolute mean flows, the integrals in (5.10) are 
zero if p ^ u, since the radial shape functions are orthogonal. This is not true, in general, 
for the acoustic perturbations of arbitrary mean flows. However, if the modal integrals do 
not vanish for p / v, then ((ab)) in (5.10) will be made up of terms that vary harmonically 
with £. In this case, we can determine the axial mean value of ((56)) as 

<«»»> = tfAi'Re ( £ K r \ 2 f° r^A’I’L,)' rdr\ . (5.11) 


Sound Pressure Level 

If we set a and 6 both equal to the unsteady pressure, p, in the foregoing equations, then 
56 = p 2 is the square of the pressure, (56) = (p 2 ) is the mean-square pressure, ((56)) = ((p 2 )) 
is the area-averaged, mean-square pressure, and 

«(ai») = «(? 2 )» - , 2 1 , , E l<wl 3 r^ rdr (5 ' 12) 

~ Vh) m,/i PTOp Jt " 
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is the axial mean value of the area- averaged, mean-square pressure. The sound pressure 
level, SPL, of the acoustic response is based on ({{ p 2 ))) and defined as: 


SPL = 101og 10 (^^f), (5.13) 

where P Re f = 2 x 10 -5 N/m 2 = 4.17706 x 10 -7 lbf/ft 2 . 

Sound Power Level 

If we set a = p~ l (p + pV x v x ) and b = pv x + pV x , then ab is the axial component of 
the acoustic energy flux, (ab) is the acoustic power or intensity, ((ab)), is the area-averaged 
acoustic power, and 


(((ab))) = Trd^Rel 



^ m,/z Pr °P 



(5.14) 


is the axial mean value of the area-averaged acoustic power. The sound power level, PWL, 
of the acoustic response is defined by 


PWL = 101og 10 


( A D {({Sb))) \ 
V PWm ) 


(5.15) 


where PW Ref = 1(T 12 watts = 0.737572 x 1(T 12 ft-lb f /sec. 

5.3 Wake Excitation Model 

The wake-excitation model, currently used with LINFLUX, is based on classical, thin, 
viscous-layer concepts, and on the foregoing far-field analysis. Thus, it is assumed that the 
flow within each wake is parallel and changes in the pressure and density across each wake are 
negligible. In addition, a quasi-3D approximation is used to determine the axial behavior, 
at inlet, of the “convected” wake perturbation velocity. 

Consider a reference rotating or stationary blade row, cf. Figure 2, which receives the 
flow exiting from an adjacent upstream row; i.e., a stator or a rotor, respectively. The 
upstream and reference blade rows rotate at constant angular velocity Q,' and Q. Therefore, 
the upstream blade row rotates at angular velocity O' — Q = CIr = relative to the 

reference array, where is the absolute angular velocity of the rotor. Because the blade row&3 
move relative to each other, any downstream, circumferential, flow nonuniformity produced 
by the upstream blade row; in particular, one due to the viscous blade wakes, will excite an 
unsteady flow through the reference row. Although the wake excitation is due to viscous 
effects, it is modeled here as an inviscid velocity perturbation at the inlet to the responding 
or reference blade row. 

We will analyze the flow at the inlet to the reference blade row in terms of cylindrical 
coordinates (r, 9, f) and (r, O' = 6 - Cl R t, £) attached to the reference and the upstream or 
excitation blade rows, respectively. Here, 9 and S' measure angular distance in the e ( xe r = e e 
or clockwise direction, when viewed from upstream, and the radial unit vector e r points 
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from hub to tip. We assume that in the absence of the reference blade row, the flow “far” 
downstream of the excitation blade row is steady in the r, 9', £ coordinate frame, with relative 
velocity, V'(r, 0',£), density, p, and the pressure P, given by 

V'M',0 = V'(r) + v R (r,0', £) + ••• 

= W'{ r ) + v R (r, 6', £)]©t(p) + ■ • ■ =W(r,^,Oer(0. (5.16) 

p(r,0',£) = p(r) + ... and P(r, 9 ', £) = P{r) + . . . . 

Here V r , p and P are circumferentially-averaged flow variables, v R (r, 9, £) is the first-order 
velocity perturbation due to the wakes from the upstream blade row, the density and pressure 
fluctuations associated with this velocity perturbation are assumed to be negligible, and e T (r ) 
is a unit vector tangent to the cylinder r — constant and pointing in the far-downstream 
flow direction. 

Thus, far downstream of the excitation blade row, the radial component of the total fluid 
velocity is negligible and the circumferentially-averaged fluid velocity, V', depends only on 
radius. The circumferentially-averaged flow is assumed to be a solution of the conservation 
equations; i.e., the nonlinear Euler equations, for steady inviscid flow. However, the wake 
perturbation will be modeled as an approximate solution of the linearized Euler equations. 
Note that, if the upstream blade row is a stator, ft' = 0, V' is the absolute steady velocity, 

12 = —ft R e^, and V is the velocity relative to the reference rotor blade row. If the upstream 
blade row is a rotor, ft' = ft R e^, V' is the velocity relative to the rotor, ft = 0, and V is the 
absolute fluid velocity. In both situations, S2 7 — 12 = fluG^, where Q,r is the absolute angular 
speed of the rotor. Also, the perturbation velocity, v R , is invariant under a transformation 
from the r, 9’,( t frame to the r, 9, £ frame. 

The flow speed, V' = V' • e T , can be represented as a Fourier series in 9'\ i.e., 

OO oo 

V'(r,9 ',£ ) = ^2 K( r -Oexp(-mI\T'0') = V^r) + ^Re{u R (r,£) exp(-mW)}, (5.17) 

7l=— CO 71—1 

where N' > 0 is the number of blades in the upstream blade row and 

l\Tf p6'+2n/N r 

V;(r, 0 = 7T ^'( r > 0 '> 0 exp n = 0, ±1, ±2, ... . (5.18) 

2 tt Jp 

Here, Vq = V' is the circumferentially-averaged flow speed far downstream of the excita- J 
tion blade row, and u R (r, £) = 2V^, n = 1, 2, . . . , is the complex amplitude of the nth 
circumferential harmonic of the perturbation flow speed. The velocity V” can be determined 
experimentally or from a Navier-Stokes solution for the steady flow relative to the upstream 
blade row. 

We assume that the Fourier components of the perturbation velocity, v R = u R e T , vary 
harmonically with axial distance, i.e., v^(r,^) = u R (r) exp(in/c^). Then, the fundamental 
wave number of this “convected” velocity perturbation is given by 

k = K T e T = (KfVf - r~ 1 N'Vj)e T /V' = 0 , (5.19) 
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where k t is the wave number in the relative flow direction. 

The flow exiting from the upstream or excitation blade row enters the field of the reference 
blade row. In particular, the imposed velocity at the inlet to the reference blade row is 

V = V' + xr = V' + Sl R re g + v R + . . . = V + v R + . . . . (5.20) 

Thus, in terms of the coordinates r, 0, f fixed to the reference blade row, the steady and 
unsteady perturbation velocities at inlet are given by 

V(r) = + (Vg + ti R r)eg — 4- Vgeg (5-21) 

and 

CO 

v R = ^2 Re { v n exp[in(/s € f - N'6 + N'n R t)]} , (5.22) 

71=1 

where 

v R = + Vgeg ) = {v„ /V')[V^ + ( Vg — fiK^)ee] , (5.23) 

N'Q r = co is the fundamental excitation frequency in the frame of the reference blade row, 
and = — (u — r~ l N'Vg)/V l \ is the axial wave number of the fundamental wake excitation. 
The fundamental disturbance repeats N' times around the wheel and has an interblade phase 
angle, a, of — 27r N 1 /N, where N is the number of blades in the reference blade row. 

At each radial station, the first-order velocity perturbation in (5.22) satisfies the lin- 
earized, inviscid, conservation equations for mass, axial and circumferential momenta, and 
energy. However, it does not conserve radial momentum, unless the absolute steady back- 
ground flow is uniform. Thus, for a nonuniform, absolute mean flow, it is, at best, only 
a convenient, and hopefully useful, approximate solution to the fluid-dynamic equations 
governing the linearized inviscid flow through the reference blade row. 
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6. ADP Rotor- Wake/FEGV Interactions 


We have applied the LINFLUX analysis to predict unsteady flows through the fan exit 
guide vane (FEGV) of the NASA/PW, 22 inch, advanced ducted propulsor (ADP). These 
flows are excited by circumferential nonuniformities in the flow exiting the fan rotor, which 
are due primarily to the viscous wakes extending downstream from the rotor blades. We 
will consider two ADP power settings; a low power or approach setting and a high power 
or take-off setting, and apply LINFLUX to predict the unsteady aerodynamic responses of 
the FEGV to Fourier components, at multiples of the blade passing frequency (BPF), of the 
rotor wake excitation. 

A schematic of the advanced ducted propulsor is shown in Figure 3. The fan consists of 18 
blades which rotate clockwise, when viewed from upstream, about the duct axis. The radius 
of the fan is 11 in and the nominal blade chord, measured over the outer span, is 3.5 in. The 
fan exit guide vane , cf. Figure 4, is placed at an axial distance from the rotor-tip trailing 
edge to the stator-tip leading edge of 5.3 in. It consists of 45 blades, which are twisted, 
bowed, flared at the tip, and have rounded leading and trailing edges. The FEGV resides 
in an aft duct of variable inner and outer radii. The chord, c*, and axial chord, c^, of the 
stator vanes, at midspan, are 1.666 in and 1.625 in, respectively. The chord is the reference 
length used in non-dimensionalizing the TURBO and LINFLUX variables. It is taken to be 
the linear distance between the vane leading- and trailing-edge points, at midspan. 

Measurements of the rotor exit flow, with the stator moved well aft of its design position, 
were taken at NASA Glenn Research Center [Pod97]. Circumferentially-averaged values of 
the absolute total temperature and absolute total pressure were determined at 10 radial 
stations in an axial measurement plane located 3.720 in downstream of the rotor-tip trailing 
edge. In addition, LDV measurements of the rotor exit velocity were taken at 29 radial and 
51 circumferential locations in an axial plane located 2.650 in downstream of the rotor-tip 
trailing edge. This data has been post-processsed to provide analytical inflow conditions for 
the TURBO and LINFLUX, steady and linearized unsteady flow calculations. The measured 
velocity field downstream of the fan, for the low power setting, cf. Figure 5, shows strong 
nonuniformities in the rotor exit flow due to rotor-wake and duct-wall phenomena. 

6.1 FEGV Geometry and Computational Grid 

To enable the application of the TURBO and LINFLUX analyses, we have introduced 
several geometric modifications to the actual FEGV configuration. In particular, a wedge-^ 
shaped, trailing-edge sections, approximately 0.2 in long, have been added to the original 
vanes, so that, the analytical vanes close in sharp trailing edges. This modification is needed 
to eliminate time-dependent, trailing-edge separations from the TURBO predictions for the 
nonlinear background flow; thereby allowing the TURBO calculations to converge to a steady 
solution. In addition, we have modified the aft-duct geometry, to provide a duct with inlet 
and exit sections of constant inner and outer radii, in accordance with the assumptions used 
in developing the steady and unsteady far-field analyses. 

The computational domain prescribed for the TURBO and LINFLUX, FEGV calcula- 
tions extends axially 1.593 in upstream and 3.407 in downstream from a blade leading edge 
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at midspan. The constant radii duct sections at inlet and exit extend almost to the blade 
row and have inner and outer radii of r* H = 5.810 in and r* D = 11.410 in at inlet; and of 
r* — 6.426 and r* D — 11.288 in at exit. In terms of the dimensionless coordinates used in the 
TURBO and LINFLUX calculations; i.e. £ = £*/e* and r = r*jc\ the key axial coordinates 
are £ = -0.956 at the computational inlet boundary, £ = 0 at the midspan leading-edge 
point, £ = 0.975 at the midspan blunt trailing-edge point, £ = 1.085 at the midspan sharp 
trailing-edge point, and £ = 2.045 at the computational exit boundary. The inner and outer 
radii, r = r*/c*, are 3.487 and 6.849 for the inlet duct section, and r H = 3.857 and r D = 6.776 
for the exit section. The axial plane at which the total temperature and total pressures were 
measured occurs at £ = -0.948; that at which the velocity measurements were taken, at 
£ = £e = —1-591. 

The TIGER grid generator [SS91] has been used to define a three-dimensional H-gnd, for 
the FEGV steady and unsteady flow calculations. The numerical results, reported herein, 
were determined on an H-grid, consisting of 173 axial, 25 radial, and 41 tangential surfaces 
(165,120 cells). This grid extends axially from £ = -0.956 to £ = 2.045 with 67 axial grid 
surfaces positioned upstream of the blade row, 75 intersecting the modified blade surfaces 
and 31 positioned downstream of the blade row. Axial grid surfaces are clustered near the 
blade leading and trailing edges; circumferential surfaces, near the blade suction and pressure 
surfaces; and radial surfaces, near the blade tip. 

Four radial stations, which we will refer to as reference stations, of the H-grid used for 
the FEGV calculations are shown in Figure 6. The reference radial stations occur near 
the tip, near seventy five percent span, near midspan, and near the hub. The LINFLUX 
near-field, finite-volume solutions, computed on this grid, are coupled to far-field acoustic 
eigensolutions, which have been determined on a radial grid consisting of 36 points clustered 
near the hub and duct walls. For the unsteady flows under consideration, the unsteady wake 
excitations are prescribed at inlet, and convected and/or nearly convected disturbances are 
numerically convected through the computational outflow boundary. 

6.2 Rotor- Exit /Stator-Inlet Conditions 

As in §5, we assume that the flow downstream of the fan can be represented as the sum 
of a circumferentially-averaged steady relative flow, in which the flow variables vary only 
with radius, and a first-order perturbation, in which the velocity is aligned with the relative 
mean flow velocity and the density and pressure are negligible. Moreover, in an attempt to 
prescribe an excitation that is an approximate solution of the unsteady equations, we require 
that the perturbation velocity satisfy a two-dimensional convection equation at each radiala 
station. Thus, the absolute steady and unsteady velocities, cf. equations (5.21) and (5.22), 
at the stator inlet have the form 

V abs = U rel ex + dreg = -f (V^ rel + £lr)eg (6-1) 

and 

v R = v R er = ^ Re exp[m(/c ? £ - AT B 0 abs -h 7V B Qt)]} , (6.2) 

71=1 

respectively. The superscripts abs and rel refer to quantities measured in the stator and 
rotor frames, respectively, u; abs = N n Vl, k ( = -(u; abs - N B V d &bs )/V^ and <r = -2nN B /N v 
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are the (reduced) frequency, axial wave number and interblade phase angle, respectively, of 
the fundamental or 1 BPF wake excitation, iV B and Ny are the numbers of rotor blades and 
stator vanes, respectively, and U is the angular speed of the rotor. The flow velocities V rel 
and v„, n = 1,2, . . . , can be determined by a Fourier analysis of the rotor exit velocity, 
cf. (5.17) and (5.18). The state vector, u c for this quasi-3D, rotor- wake excitation; i.e., 
(uc) T = p[ 0, 0, vf, uf, l/ e abs {;^ + describes a first-order perturbation in which mass, 

axial and circumferential momenta, and energy are conserved, but radial momentum is not 
conserved, unless V abs is uniform. 

Analytic Inflow Conditions 

The inflow conditions for the FEGV, steady and unsteady calculations; i.e., the radial 
distributions of the circumferentially-averaged, absolute, total temperature, T£ bs , total pres- 
sure, P| bs , and tangential flow angle, fi abs = tan -1 and the complex amplitudes, 

(r), of the wake perturbation velocity, are determined by post-processing the NASA data 
for the rotor exit flow. Post-processing is required to avoid inlet conditions that lead to 
separation in the TURBO solution; remove axial vorticity phenomena at the duct walls and 
provide a more or less classical wake excitation that can be accommodated by LINFLUX; 
and to account for the diffusion of the viscous wakes as they are convected downstream. 

The analytic inlet, total temperature and total pressure distributions are obtained by 
fitting cubic polynomials to the NASA data; the analytic flow angle distribution, by fitting 
a linear curve for r H < r <r M and continuing this into a quadratic curve for r M < r < r D , 
where tm is the approximate radial location at which fi abs reaches a minimum. The fitting 
coefficients are adjusted to avoid inflow angles that lead to separation in the TURBO solution. 

The analytic wake excitation has been derived from the NASA velocity data in the fol- 
lowing manner. First, the measured wake velocity perturbation is determined by subtracting 
the circumferentially-averaged values of the measured, relative velocity from the actual rel- 
ative velocity. Since the radial velocity was not measured, and is assumed to be zero, this 
yields a velocity perturbation with components, that are parallel, vt, and normal, 5 tv, to the 
circumferentially- averaged, relative velocity at the velocity measurement plane, £ = £e- The 
normal component of the wake perturbation velocity is generally found to be small. Thus, 
in keeping with the assumptions used in developing the LINFLUX wake-excitation model, 
only the parallel component of the measured perturbation velocity is used to construct the 
analytic wake excitation velocity. 

In current wake-excitation models [MG84, TE99, TopOO], the wakes are assumed to be 
identical from blade to blade. At each radial station, the perturbation velocity, vt{N), in^ 
a single wake, is assumed to be symmetric relative to a wake centerline at N — 0. The 
velocity distribution, is described in terms of the wake-edge, 5r,e, and centerline, 

v T , min, velocities and a shape function, f{N). In the present study, a symmetric, hyperbolic 
secant, velocity profile has been fitted to the NASA data for fir(r, 6 , to provide an analytic 
description of the wake excitation velocity at each radial station. In addition, the radial and 
circumferential distributions of the fitted wake velocity are modified near the duct walls to 
remove splitter and end wall effects (cf. Figure 5), and provide a more or less classical wake 
excitation from hub to tip. The fitted wake perturbation velocity is assumed to act in the 
direction of the analytic, rather than the measured, relative mean velocity. 
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Finally, the diffusion of the viscous wakes, as they approach the FEGV must be taken 
into account. Based on the fitted velocity profiles, the wake velocity defect, v T>e - v Tiinin , 
and half-width at the measurement plane are determined as functions of radius. Empirical 
relations [MG84] are then applied to describe the effects of viscous diffusion on the behavior of 
the wake velocity defect and half-width, and hence, on that of the wake velocity distribution, 
with increasing streamwise distance downstream of the rotor. The viscous wake velocity 
profiles at some reference axial plane, taken here to be the FEGV, midspan, leading-edge 
plane, are then used to determine the Fourier coefficients, v%, of an equivalent inviscid 
wake-excitation at the FEGV, computational, inlet plane. Thus, the Fourier components 
of the inviscid wake excitation that are prescribed, at inlet, to LINFLUX, are based on 
an empirically determined, viscous, rotor-wake, velocity distribution at the FEGV, midspan, 
leading-edge plane. This viscous velocity distribution is, in turn, based on the measurements 
taken downstream of the rotor at NASA. 

Numerical Results 

At this point, we have described the analytic FEGV vane and duct geometries, the grid 
used for the steady and unsteady flow calculations, and the manner in which the analytic 
inflow conditions were determined from the NASA data. We will proceed to present numeri- 
cal results in §7 and §8 that describe the flows through the FEGV at the low and high ADP 
power settings. AWAKEN results will be given for the circumferentially- averaged inflow 
conditions, the rotor-wake excitation velocity at the FEGV leading edge, and the Fourier 
components of this velocity. TURBO predictions will be given for the steady pressure field 
and the steady vane-surface pressure distributions. Finally, LINFLUX results will be given 
for the imaginary components of the unsteady velocity and pressure fields; the real and imag- 
inary unsteady vane-surface pressure distributions; the axial eigenvalues of far-field acoustic 
disturbances; the radial shapes of the propagating acoustic response modes; and the modal 
and total, sound pressure and sound power levels of the acoustic responses at inlet and exit. 
The predicted sound power levels at exit will be compared with measured data. 

The component codes of the LINFLUX Prediction System use different forms of the flow 
variables. As a result, the rotor-exit/FEGV inflow conditions are reported below, in terms 
of the dimensional flow variables used in AWAKEN. The steady and unsteady flows through 
the FEGV are, for the most part, described in terms of the non-dimensional flow variables, 
discussed in §2 of this report, and used in the LINFLUX output routines. Steady inflow and 
outflow results are also reported in terms of the non-dimensional variables used in TURBO. 

For the low power setting, unsteady response predictions have been determined for rotor- 
wake excitations at 1, 2, and 3 times the blade passing frequency. However, the results for 3 
the 3 BPF excitation are inaccurate. They are included here for the sake of completeness 
and for possible future reference. Unsteady results at the high power setting are given for 
1 and 2 BPF wake excitations. In this case, a converged LINFLUX solution could not be 
determined for the 3 BPF excitation. The TURBO and LINFLUX calculations have been 
performed on the 173 x 25 x 41 near-field mesh, shown in Figure 6. The LINFLUX far-field 
analysis was performed on a 36 point radial mesh. The TURBO calculations required several 
days of CPU time on an IBM 3CT workstation to converge to an attached, low-loss, steady 
solution. Typically, the LINFLUX calculations required 20 CPU hours on the same machine 
to converge to an unsteady solution. 
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7. Numerical Results: Approach Power Setting 


At the low power setting, the fan of the ADP rotates at an angular speed of 5,425 rpm and 
the mass flow rate through the FEGV is 45.2 lb m /sec. Radial distributions of the measured 
and analytic, absolute, total temperature, T| bs , total pressure, P£ hs , and tangential flow 
angle, Q abs = tan -1 (V^/V^), at the rotor exit (stator inlet) are shown in Figure 7. The 
analytic distributions are used as inlet conditions for the TURBO calculation of the steady 
background flow through the FEGV. Note that, in the tip region, the analytic inflow angles 
have been set deliberately at values less than those measured, so that the TURBO analysis 
will yield an attached steady-flow solution. The analytic exit-plane static pressure at the 
hub is determined by adjusting the prescribed value of this pressure, by trial and error, until 
the analytical and experimental mass-flow rates arc the same. 

Although the inflow values of T| bs and P| bs , prescribed for the TURBO calculation, are 
in excellent agreement with the corresponding experimental data, and the values of ft abs are 
in good agreement, there are differences between the calculated and the measured mean flow 
velocities. The circumferentially-averaged, measured, relative and absolute velocities at the 
rotor exit (stator inlet) are shown in Figure 8, along with the values determined from the 
TURBO steady flow solution. The differences between the measured and calculated steady 
inflow velocities are due, in part, to the differences between the measured and prescribed 
values of fl abs . However, the attempt to model the measured viscous flow by an analytical 
inviscid flow at the same mass flow rate is probably a more important contributor. 

Circumferential distributions, over 3 rotor passages, of the analytic, rotor-wake, velocity 
perturbation, v R = v R er, at the midspan leading-edge plane of the stator, as determined 
from the NASA data and the post-processing steps outlined in §6, are shown in Figure 9. 
In this figure, 0 G is the angular gap of the rotor and, for clarity, the velocity profiles at each 
radius are offset by rjr D x 10 3 . The curves in Figure 9 indicate only the behavior of the 
magnitude of the tangential perturbation velocity, but there are also strong radial changes 
in the direction, e T , of this velocity, as indicated by the solid curve for Q rel in Figure 8. The 
circumferential location at which the minimum wake velocity occurs, also varies significantly 
along the span. The velocity profiles in Figure 9 indicate that the wake velocity defect 
and half-width are relatively large near the hub and tip, and much smaller over the blade 
midspan region. Some of this behavior could be due to endwall and splitter effects that are 
still present in the post-processed wake velocity data. 

Radial distributions of the amplitude, |v*|, and phase, arg(v R ), of the first five Fourier 
components of the analytic wake excitation at the stator leading edge are shown in Figure 10.^ 
The amplitudes of the Fourier excitations decrease with increasing circumferential mode 
number, n. Also, the amplitude of the first harmonic excitation shows a strong radial 
dependence, varying from approximately 25 ft /sec at the hub, to 10 ft/sec at midspan, to 
27 ft/sec at the tip. The higher harmonics show much smaller variations in amplitude. In 
particular, the amplitude of the 2 BPF excitation varies from 7 to 10 ft /sec over the span; 
that of the 3 BPF excitation, from 3 to 5 ft/sec. Large radial variations in phase also occur; 
e.g., a phase angle variation of approximately 160 deg for the first harmonic, 360 deg for the 
second harmonic, and 520 deg for the third harmonic. 
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7.1 Steady Background Flow 

The steady flow through the FEGV, for the low power setting, has been determined by 
applying the TURBO 3D Euler analysis on the computational grid illustrated in Figure 6. 
The radial distributions of absolute total temperature, absolute total pressure and absolute 
flow angle, indicated by the solid curves in Figure 7, are prescribed at the computational 
inlet plane, located at f = -0.956, and the static pressure at the hub is specified at 
exit, £ = £ + = 2.045, such that the calculated mass-flow rate matched that measured in 
the NASA Glenn experiment. The computed results indicate that the isentropic, absolute, 
Mach number at inlet increases from 0.240 at the hub (r = 0.509r£>) to a maximum of 0.346 
at t = 0.827ro, and decreases from 0.346 to 0.309 at the tip. At the exit, this Mach number 
increases from 0.151 at the hub (r = 0M9r D ) to a maximum of 0.318 at r = 0.866^, and 
decreases from 0.318 to 0.294 at the tip. 

Selected results from the TURBO calculation are given in Figure 8 and in Figures 11 
through 13. The steady pressure field is shown, at the four reference spanwise stations, in 
Figure 11. The pressure is nearly constant at the computational inlet and exit boundaries, 
and the steady pressure gradients, generated by the flow through the FEGV , die out within 
short axial distances from the blade row. Because of steady blade loading, there is a small 
pressure rise across the FEGV. The steady pressure distributions over the suction and pres- 
sure surfaces of a modified vane, i.e., the original vane with a wedge-shaped trailing-edge 
section added, are shown in Figure 12. The pressure distributions indicate an increase in 
blade loading from the hub to approximately 75% span. There are flow overspeeds, just 
downstream of the vane leading edge, at the J =16 and J =22 stations, leading to peak 
Mach numbers of 0.499 and 0.435, respectively, on suction surface of each vane. Aft of 
the leading-edge region, peak isentropic Mach numbers of approximately 0.359, 0.469, 0.469 
and 0.440 occur in the vicinity of £ = 0.25^ at the J = 4, 10, 16 and 22 radial stations, 
respectively. 

Radial distributions of the steady density, pressure, and velocity at the computational 
inlet and exit planes of the FEGV are shown in Figure 13. To separate the curves in Figure 13, 
the results are reported in terms of the non-dimensional flow variables, used in the TURBO 
code. Thus, the density is scaled by the reference density, p^ f = 0.0764 lb m /ft 3 ; pressure, 
by the reference pressure, P£ ef = 2, 116.2 lb f /ft 2 ; and velocity by the reference sound speed 
^Ref — 1,116.8 ft /sec. The reference density and sound speed are determined from the 
reference pressure and temperature, T£ ef = 519 °R, using the perfect gas relations. 

The results in Figure 13 indicate that the steady pressure and density are nearly constant 
over the inlet and exit planes, and there is a slight rise in these state variables across the_j 
blade row. The axial velocities at inlet and exit have different radial distributions, with a 
greater radial variation occurring in the axial velocity at exit. The circumferential veloc- 
ity distributions indicate that the steady flow at inlet has appreciable swirl, with the swirl 
velocity varying from 0.163A^ ef at the hub to 0.208A^ ef at r = 0.74 6r D . The FEGV essen- 
tially removes this swirl, leading to a nearly axial exit flow. Finally, the inlet flow has, by 
prescription, zero radial velocity, but the exit flow has a small radial velocity component. 
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7.2 Unsteady Response Predictions 


We have determined the unsteady aerodynamic responses of the FEGV, for the ADP 
operating at the low power setting, to rotor wake excitations at 1, 2, and 3 BPF. These 
excitations occur at reduced frequencies, u = nN B Cl, of 3.703, 7.406 and 11.109, respectively, 
and interblade phase angles, a = —2n nN B fNy, of -144 deg, -288 deg and -432 deg. They 
repeat nN B times around the wheel and travel toward the FEGV and in the direction of 
rotor rotation. In addition to occurring at moderate to high reduced frequencies, the wake 
excitations are highly dependent on radius. These features lead to rather complicated three- 
dimensional aerodynamic responses. The FEGV unsteady response predictions for the low 
ADP power setting are given in Figures 14 through 28. 

The LINFLUX far-field eigenanalysis determines the axial and radial behaviors of the 
modal acoustic disturbances far upstream and far downstream of the FEGV. These be- 
haviors depend upon the underlying steady background flows at inlet and exit, the ex- 
citation frequency and interblade angle, and the circumferential angular wave numbers, 
77 i = —nN B + mNy , of the modal unsteady responses; but, not on the type, amplitude, 
or radial behavior of the unsteady excitation. In the present effort, the far-field eigenanaly- 
sis has been applied to determine the axial eigenvalues and radial shapes of the lowest, 
radial-order, acoustic modes in the m = -2, -1,0, 1,2 circumferential modes, for the 1 BPF 
excitation, and in the m = — 3, ..., 3 circumferential modes, for the 2 and 3 BPF excitations. 
Although Some modes can be missed or filtered out by the far-field analysis, usually, such 
modes are highly attenuated and have little impact on the overall unsteady solution. 

We indicate the axial eigenvalues of both acoustic excitations and acoustic responses 
in our eigenvalue plots. At inlet, the upstream propagating or decaying disturbances are 
acoustic responses, and the downstream propagating or decaying disturbances are acoustic 
excitations. At exit, the reverse is true. For the present application to wake/blade-row 
interactions, the amplitudes of the modal acoustic excitations are set to zero. The amplitudes 
of the modal acoustic responses are determined by the coupling of the LINFLUX near- and 
far-field solutions. 


1 BPF Wake Excitation 

The axial eigenvalues of the modal acoustic disturbances in the far-field of the FEGV for 
the 1 BPF excitation are shown in Figure 14, where we have indicated the circumferential, m, 
and radial, /*, mode orders of the the least-damped acoustic response disturbances. For the 1 
BPF excitation, all acoustic response modes attenuate. Moreover, because of the swirl in the ^ 
steady background flow at inlet, there are differences between the axial wave numbers of the 
attenuating disturbances that occur in different circumferential modes. At 1 BPF, the least- 
damped response modes at inlet and exit occur at m, n = 0, 0. These repeat eighteen times 
around the wheel; travel in the direction of rotor rotation (fh = —nN B + mNy = —18); 
have no radial nodes (/* = 0); and attenuate with increasing distance from the FEGV. 
The axial eigenvalues of the upstream and downstream traveling 0,0 acoustic responses are 

Xo.o = 2-924 + 0.209* and Xo,o = -2.786 + 0.443*. 

' The FEGV unsteady velocity and pressure fields resulting from the 1 BPF wake excita- 
tion are shown in Figures 15 and 16, respectively. Because plots of the real and imaginary 
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components of the unsteady field variables indicate similar physical phenomena, only the 
imaginary components of the unsteady velocity, 7m{v}, and pressure, Im{p}, will be de- 
picted in this report. The unsteady velocity vectors, in Figure 15, have been obtained by 
multiplying the predicted velocities by a scale factor of 12. The velocities at the computa- 
tional inlet plane are , in general, due to the prescribed wake excitation and the acoustic 
response of the FEGV. However, for the 1 BPF excitation, the pressure response diminishes 
strongly with increasing axial distance from the blade row. Thus, the velocity at the in- 
let plane in Figure 15 is almost entirely due to the wake excitation. The 1 BPF unsteady 
velocity field shows a strong response at the leading-edge of each blade, and a change in 
the unsteady flow into one dominated by counter-rotating vortical eddies, in which the flow 
speeds are low, as it proceeds through the stator blade row. 

As noted, the excitation at 1 BPF produces attenuating acoustic responses upstream 
and downstream of the FEGV. This is clearly indicated by the pressure contours plotted in 
Figure 16, where the small pressure variations at the computational inlet and exit boundaries 
are due to acoustic responses that die out with increasing distance from the blade row. The 
pressure amplitudes of the least-damped 0,0 acoustic responses, at the computational inlet 
and exit boundaries, are 0.00140 and 0.000689, respectively. The average pressure at the 
FEGV inlet is 6.101. Thus, the unsteady pressure amplitude of 0.0014 is only 0.023% of the 
average inlet pressure. The next strongest mode, in terms of amplitude, at inlet is the 1,0 
mode, with amplitude 0.00106. The next strongest at exit is the 0,2 mode, which has an 
amplitude of 0.000181. However, the inlet 1,0 and the exit 0,2 modal responses are more 
highly damped than the inlet and exit, 0,1, response modes. 

The far-field pressures are the response quantities of primary interest in blade-row aeroa- 
coustic response studies. However, we also present surface-pressure response predictions 
to assist in demonstrating the capabilities and accuracy of the LINFLUX analysis and to 
provide additional response information, that can be tested experimentally. The unsteady 
surface pressure distributions acting on the reference or zeroth FEGV blade for the 1 BPF 
wake excitation are shown, at the J = 4, 10, 16, and 22 radial stations, in Figure 17. Here, 
the surface pressure responses are given for the modified FEGV blades, which extend from 
the blade leading edge near £ = 0 to the sharp trailing edge near £ = Lllc^. The trailing 
edge of the original (blunt-edged) blade occurs near £ = c^. Since the amplitude, direction 
and phase of the 1 BPF wake-velocity vary significantly from hub to tip, and this excitation 
occurs at a moderate reduced frequency, it is difficult to discern any simple radial trends 
from the surface pressure responses. 

2 BPF Wake Excitation 

As indicated by the axial eigenvalues, plotted in Figure 18, the 2 BPF wake excitation 
generates propagating acoustic responses upstream and downstream of the blade row. In 
particular, three response modes, the 1,0, 1,1, and 1,2 modes, propagate upstream and two, 
the 1,0 and 1,1 modes, propagate downstream. The axial wave numbers of the 1,0 , 1,1 and 
1,2 upstream propagating modes, which travel against the incoming steady flow, are 3.387, 
2.862 and 2.049, respectively. Those of the 1,0 and 1,1 downstream propagating modes, 
which travel with the exiting mean flow, are -1.378 and -0.760. The least-damped upstream, 
1,3, and downstream, 1,2, attenuating acoustic responses, at 2 BPF, have axial eigenvalues 
XY ,3 = 1-631 + 0.765* and xt ,2 = -1-H7 + 0.702*. 
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The imaginary components of the unsteady velocity and pressure resulting from the 2 
BPF excitation are shown in Figures 19 and 20, respectively. To obtain the results in Fig- 
ure 19, a scale factor of 30 was applied to the predicted velocities. The velocity distribution 
along the computational inflow boundary consists primarily of the prescribed wake-excitation 
velocity, but there is also a contribution from the velocity associated with the propagating 
acoustic response. Again, there are high velocity responses at the leading edge of each vane 
and a counter-rotating, small-eddy, flow structure extending downstream from the vane 
leading-edge plane. 

The pressure contours in Figure 20 show pressure variations in the vicinities of the compu- 
tational inflow and outflow boundaries that are due to propagating acoustic response waves. 
The persistence of such waves is clearly indicated upstream, at the 94% spanwise station, 
and downstream, at the four reference spanwise stations. A spurious pressure behavior also 
occurs near the inlet boundary, which is probably due to the errors associated with imposing 
a wake excitation that is only an approximate solution of the governing unsteady aerody- 
namic equations. The predictions indicate complicated acoustic responses, both upstream 
and downstream of the blade row. In addition to the complex radial behavior and the rel- 
atively high reduced frequency of the wake excitation, the modal acoustic responses also 
introduce strong radial gradients. 

For the 2 BPF excitation, the amplitudes of the upstream propagating 1,0, 1,1 and 1,2 
acoustic response modes at inlet are 0.000775, 0.000177, and 0.00221, and the least-damped 
1,3 mode has an amplitude of 0.000830. Thus, the 1,2 mode occurs at the highest amplitude, 
and this is 0.036% of the average inlet pressure. At exit, the propagating 1,0 and 1,1 modes 
have amplitudes of 0.00193 and 0.00157, and the least-damped (1,2) mode has an amplitude 
of 0.000921. There is no clearly dominant acoustic mode at the computational inlet or exit 
boundaries. In addition, the modes that contribute to the pressure responses at inlet and exit 
have strong radial variations. Such features make it difficult to provide a simple or convenient 
interpretation of the far-field acoustic responses, based on contours plotted in Figure 20. An 
additional complication is introduced by the spurious pressure behavior at inlet. Note that 
the pressure contour levels for the 2 BPF excitation, in Figure 20, are much lower than 
those for the 1 BPF excitation, in Figure 16. The 1 BPF excitation produces a stronger 
unsteady pressure response near the blades, but no sound in the far field. Therefore, the 
2 BPF excitation is the more important one aero-acoustically, because it produces pressure 
responses that persist far from the blade row. 

The real and imaginary unsteady surface-pressure responses to the 2 BPF wake excitation 
are shown in Figure 21. Note that the (ordinate) scale used to plot the surface- pressure 
response for the 2 BPF excitation is roughly one-fifth of that used for the 1 BPF excitation .-3 
The real and imaginary surface pressure responses to the highly three-dimensional 2 BPF 
wake excitation show no discernable simple trend, but, instead, a relatively complicated 
radial behavior. 

The radial shapes of the propagating acoustic response modes, produced by the 2 BPF 
excitation, are shown in Figure 22. The propagating response disturbances all occur in the 
m = 1 circumferential mode, i.e., at m = —2 Ng + Ny = 9 or a = 2irm/Ny = 72 deg. Hence, 
they repeat nine times around the stator blade row, and travel circumferentially counter to 
the direction of rotor rotation. The propagating modes, including those at // = 0, all show 
strong spanwise pressure variations. 


NASA/CR— 2001-210713 


29 



We have determined the sound pressure and sound power levels of the 2 BPF acoustic 
responses. The predicted sound pressure levels, in dB, for the 1,0, 1,1 and 1,2 acoustic 
response modes at inlet are 110.4, 97.0 and 117.0, respectively. Those for the 1,0 and 
1,2 response modes at exit are 117.6 and 115.7. The total sound pressure levels at inlet 
and exit are 117.9 and 119.8, respectively. The predicted modal and total sound power 
levels are given in Figure 22, along with the measured power levels (in parentheses), at 
exit [Hei99]. Considering the various approximations made at this stage of the development 
of the LINFLUX Prediction System, particularly those associated with modeling the wake 
excitation, the level of agreement between the predicted and measured data is encouraging. 

It is also interesting to note that pressure responses at relatively small amplitudes; i.e., at 
amplitudes that are much less than 0.1% of the average inlet pressure, result in sound power 
and sound pressure levels on the order of 115 dB. 

3 BPF Wake Excitation 

As mentioned previously, the 3 BPF response predictiions are not accurate. They are 
included here for future reference, and because they illustrate important capabilities and 
limitations of the LINFLUX code. The 3 BPF wake excitation produces an unsteady aero- 
dynamic response in which, as indicated in Figure 23, four acoustic response modes, the 1,0, 
1,1, 1,2 and 1,3 modes, propagate upstream and downstream of the FEGV. The axial wave 
numbers of the 1,0, 1,1, 1,2 and 1,3 upstream propagating modes are 4.441, 3.846, 3.190 and 
2.009, respectively; those of the corresponding downstream propagating modes are -2.574, 
-2.372, -1.782 and -0.369. The least-damped, upstream and downstream, attenuating re- 
sponses at 3 BPF have axial eigenvalues of X 14 = 2.252 + 0.924i and xtp = -4.347 + 1.276L 
A damped 1,4 mode was not captured by the downstream far-field analysis. 

The unsteady velocity field associated with the 3 BPF rotor-wake/FEGV interaction is 
shown in Figure 24, where the predicted velocities have been multiplied by a scale factor of 
75. The velocities at the inlet plane are due primarily to the prescribed wake excitation, but 
departures from the prescribed velocities, due to the response of the FEGV, are evident from 
the from the behavior of the velocity vectors in the vicinity of the inlet plane. Downstream of 
the vane leading edge the unsteady flow speeds become quite small, and the velocity vectors 
are indicative of an unsteady flow that is dominated by counter-rotating, vortical eddies. 

The unsteady pressure field generated by the 3 BPF wake/FEGV interaction is depicted 
in Figure 25. The contours in Figure 25 indicate a very complicated acoustic response 
upstream of the FEGV, and a spurious response near the computational inlet boundary. 
There is, however, a well-defined propagating acoustic response downstream of the blades 
row. The 3 BPF wake excitation varies strongly with radius. It also occurs at a high 
frequency and at small amplitude. Thus, prescribing this excitation, in terms of a quasi-3D 
wake model, probably leads to errors that are significant relative to the true aerodynamic 
response of the FEGV. Such errors are especially severe near the inflow boundary. 

The far-field acoustic eigensolutions, i.e., the axial eigenvalues and radial mode shapes, 
depend upon the frequency and circumferential wave number of the unsteady excitation, but 
not on the amplitude or radial content of this excitation. The latter do, however, impact the 
amplitudes of the modal acoustic responses. The LINFLUX predictions indicate that the 
amplitudes of the upstream propagating 1,0, 1,1, 1,2 and 1,3 modal acoustic responses at 
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inlet are 0.000505, 0.000278, and 0.000369, and 0.000712, respectively. The strongest damped 
mode; i.e., the 0,2 mode, has an amplitude of 0.000151. Surprizingly high amplitudes are 
predicted for the more heavily damped modes, which, again, indicates error in the unsteady 
solution. At exit, the propagating 1,0, 1,1, 1,2 and 1,3 modes have amplitudes of 0.000207, 
0.0000607, 0.000370 and 0.000122, and the damped modes are quite weak. There are no 
clearly dominant acoustic modes at inlet or exit. The inlet, 1,3 response occurs at the 
highest amplitude, and this is only 0.012% of the average inlet pressure. Thus, the modes 
that contribute to the pressure responses at inlet and exit occur at small amplitude, have 
strong radial variations, and reach their maximum amplitudes at different radial locations. 

The unsteady surface-pressure responses of the FEGV to the 3 BPF wake excitation 
are shown in Figure 26. The real and imaginary surface pressure responses to this highly 
three-dimensional, excitation show relatively complex radial and axial variations. Note that 
the pressure (ordinate) scale in Figure 26 is approximately 20% of that used for the 2 BPF 
surface-pressure response (Figure 21), and only 4% of that used for the 1 BPF response 
(Figure 17). 

The radial shapes of the propagating acoustic response modes, at inlet and exit, for the 
3 BPF excitation are shown in Figures 27 and 28, respectively. The propagating modes at 3 
BPF occur at an angular wave number of m = -3 N B + N v = -9. Hence, as in the 2 BPF 
case, the model acoustic responses repeat nine times around the wheel; but, the responses 
at 3 BPF travel in, rather than counter to, the direction of rotor rotation. This difference 
in the propagation behaviors of the 2 and 3 BPF acoustic responses is also indicated by the 
orientations of the downstream pressure contours , in Figures 20 and 25. 

The predicted sound pressure levels, in dB, for the inlet 1,0, 1,1, 1,2 and 1,3 response 
modes are 105.5, 100.5, 102.6, and 107.0. Those for the 1,0, 1,1, 1,2 and 1,3 response modes 
modes at exit are 99.1, 86.4, 101.5 and 92.1, respectively. The total sound pressure levels 
at inlet and exit are 110.6 and 103.0, respectively. Modal and total, sound power level 
predictions, for the acoustic responses at inlet and exit, are given in Figures 27 and 28. For 
the 3 BPF excitation, the predicted sound power levels at exit differ substantially from the 
measured values [Hei99], which are given in parentheses in Figure 28. It should be noted 
that a downstream propagating 1,3 response at exit was not detected in the experiment. 
It is also important to mention that the validity of the experimental results for the 3 BPF 
excitation is in question. 


-3 
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8. Numerical Results: Take-off Power Setting 


We proceed to consider the steady and unsteady flows through the FEGV for the ADP 
operating at the high power setting. At this setting, the fan rotates at 8,750 rpm and the 
mass flow rate through the FEGV is 72.8 lb m /sec. The radial distributions of the measured 
and analytic, absolute, total temperature, T^ hs , total pressure, P£ bs , and tangential flow 
angle, S2 abs = tan -1 (V 0 abs / V^) , at the rotor exit (stator inlet) are shown in Figure 29. The 
analytic distributions are prescribed at inlet for the TURBO steady- flow calculation. We 
did not try to improve the agreement between the analytic and measured values of f2 abs , as 
this would entail a more detailed curve-fitting procedure, and some additional patience to 
ensure that TURBO converged to a low-loss, attached, steady solution. 

The inflow values of Tf hs and P£ bs , prescribed for the TURBO calculation, are higher 
than those prescribed at the low power setting, and they are in very good agreement with 
the measured data. The values of f2 abs are higher over the outer 50% span of the blading, 
and they are in fair agreement with the data. As indicated in Figure 30, there are differences, 
particularly near the hub, between the measured and the calculated mean flow velocities. 
These are due, in part, to the differences between the measured and prescribed values of 
Q abs , and to representing the measured viscous flow by an analytic inviscid flow at the same 
mass flow rate. 

Circumferential distributions, over 3 rotor passages, of the analytic, rotor-wake, velocity 
perturbation, v R = v R ex, at the midspan leading-edge plane of the stator, are shown in 
Figure 31. These distributions were determined by post-processing the NASA data for the 
high power setting. The curves in Figure 31, which are offset by r/r# x 10 3 , describe the 
behavior of the coefficient, v R , of the wake perturbation velocity. There are also radial 
variations in the direction, er, of this velocity, as indicated by the solid curve for O rel in 
Figure 30. The circumferential location at which the minimum wake velocity occurs, also 
varies significantly along the span. The velocity profiles in Figure 31 indicate that the wake 
velocity defect and half-width are relatively large near the hub and tip, and much smaller 
over the blade midspan region. Also, the wake excitation, at the high power setting, is much 
stronger than that, depicted in Figure 9, for the low power setting. 

Radial distributions of the amplitude, |u R |, and phase, arg(w R ), of the first four Fourier 
components of the analytic wake excitation velocity, at the stator leading edge, are shown 
in Figure 32. Again, the amplitudes of the Fourier excitations decrease with increasing 
circumferential mode number, n. Also, the amplitude of the first harmonic excitation shows 
a strong radial variation over the span, varying from approximately 40 ft/sec at the hub, to^ 
19 ft /sec at midspan, to 33 ft/sec at the tip. The second harmonic varies from 8 ft/sec at 
the hub to 17 ft/sec at the tip. Large phase variations also occur; e.g., a phase variation 
of approximately 210 deg for the first harmonic and approximately 420 deg for the second 
harmonic. 

8.1 Steady Background Flow 

The steady flow through the fan exit guide vane, for the high power setting, has been 
calculated by applying the TURBO analysis on the computational grid illustrated in Figure 6. 
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The radial distributions of absolute total temperature, absolute total pressure and absolute 
flow angle, indicated by the solid curves in Figure 29, are prescribed at the computational 
inlet plane, located at £ = = —0.956, and the static pressure at the hub is specified, by 

trial and error, at the exit, £ = £ + = 2.045, such that the calculated mass-flow rate matched 
that measured in the NASA Glenn experiment. The computed results indicate that the 
absolute, isentropic, Mach number at inlet increases from 0.306 at the hub (r = 0.509r D ) to 

0. 569 at r = 0.847r£>, and decreases from 0.569 to 0.517 at the tip. At exit, the isentropic 
Mach number increases from 0.144 at the hub ( r = 0.569^) to 0.519 at r = 0.898r D , and 
decreases from 0.519 to 0.488 at the tip. 

Results from the TURBO calculation are given in Figure 30 and in Figures 33 through 35. 
The steady pressure field, at the four reference radial stations, is shown in Figure 33, and has 
characteristics similar to those indicated in Figure 11 for the low power setting. In particular, 
the steady pressure is nearly constant at the computational inlet and exit boundaries, the 
pressure variations generated by the FEGV die out within short axial distances from the 
blade row, and, because of blade loading, there is a small pressure rise across the FEGV. 

Steady pressure distributions over the suction and pressure surfaces of the modified vane, 

1. e., the original vane with a wedge-shaped trailing-edge section added, are shown in Fig- 
ure 34. These indicate an increase in vane loading from hub to tip. There are flow overspeeds 
on the suction surface at J =16 and J = 22, leading to peak Mach numbers of 0.989 and 
1.130, respectively, just aft of the leading edge. Downstream of the leading-edge region, 
the peak isentropic Mach numbers at the reference radial stations are approximately 0.516, 
0.776, 0.815 and 0.767 at J = 4, 10, 16 and 22, respectively. The axial locations at which 
these peak Mach numbers occur, generally lie in the vicinity of £ = 0.2cax. 

Radial distributions of the steady density, pressure, and velocity, at inlet and exit, are 
shown in Figure 35. Here, the flow variables are scaled according to the TURBO non- 
dimensionalization. Thus, the density is scaled by the reference density, p^ ef = 0.0764 lb m /ft 3 ; 
pressure, by the reference pressure, P£ ef = 2,116.2 lb f /ft 2 ; and velocity by the reference 
sound speed ^Ref — 116.8 ft/sec. The pressures and densities at inlet and exit for the 

high power setting are slightly higher than those for the low power setting, cf. Figure 13; 
whereas, the inlet axial and circumferential velocities and the exit axial velocities are signif- 
icantly higher. 

The results in Figure 35 indicate that the steady pressure and density are nearly constant 
at inlet and exit, and there is a slight rise in pressure and density across the blade row. The 
axial velocities at inlet and exit have different radial distributions, and show strong radial 
variations, particularly in the axial velocity at exit. The circumferential velocity distributions 
indicate that the steady inlet flow has a strong swirl. The FEGV removes this swirl, leading-3 
to an nearly axial exit flow. Finally, the inlet flow has, by prescription, zero radial velocity, 
but the exit flow has a small radial velocity component. 

8.2 Unsteady Response Predictions 

We have determined the unsteady aerodynamic responses of the FEGV, for the ADP 
operating at the high power or take-off setting, to rotor wake excitations at 1 and 2 BPF. 
These excitations occur at reduced frequencies, u = nN B tl, of 3.677 and 7.355, respectively, 
and interblade phase angles, a = - 2'nnN B /N v , of -144 deg and -288 deg. Thus, the 1 and 2 
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BPF wake excitations, at the high power setting, occur at reduced frequencies similar to and 
interblade phase angles identical to those for the 1 and 2 BPF excitations at the low power 
setting. These excitations repeat nN B times around the wheel and travel axially toward the 
FEGV and circumferentially in the direction of rotor rotation. The wake excitations are 
highly dependent on radius, as indicated by the curves depicted in Figure 32. 

The LINFLUX unsteady response calculations were performed on the 173 x 25 x 41 
mesh, shown in Figure 6, and a 36 point, radial, far-held mesh. The far-held eigenanalysis 
was applied to determine the lowest-order radial acoustic modes in the m = -2, -1,0, 1,2 
circumferential modes, for the 1 BPF excitation, and in the m = -3,..., 3 circumferential 
modes, for the 2 BPF excitation. The LINFLUX results for the 1 and 2 BPF excitations at 
the ADP high-power setting are given in Figures 36 through 45. We have also constructed 
animations of these results, which are described below. 

1 BPF Wake Excitation 

The axial eigenvalues of the modal acoustic disturbances in the far held, for the 1 BPF 
excitation, are shown in Figure 36, where we have labeled the circumferential, m , and radial, 
[x, mode orders of the least-damped acoustic response disturbances. Because of swirl in 
the steady background how, the axial wave numbers of the attenuating disturbances in 
the various circumferential modes, differ. For the 1 BPF excitation, all acoustic response 
modes attenuate. The least-damped response modes are the 1,0 mode upstream, with x7o = 
2.225 + l,799i, and the 0,0 mode downstream, with x£ 0 = -2.082 + 1.420L Thus, at inlet, 
the least-damped, 1,0 acoustic response mode repeats twenty seven times around the wheel; 
travels counter to the direction of rotor rotation (m = -N B + N V = 27); has no radial nodes 
(fi = 0); and attenuates with increasing axial distance upstream from the FEGV. At exit, 
the least-damped, 0,0 mode repeats eighteen times around the wheel; travels in the direction 
of rotor rotation (m = — N B = —18); has no radial nodes {p — 0); and also attenuates with 
increasing axial distance downstream from the FEGV. 

The imaginary components of the FEGV unsteady velocity, Jm{v}, and pressure, Im{p}, 
fields, resulting from the 1 BPF wake excitation, are shown in Figures 37 and 38, respectively. 
The velocity vectors, depicted in Figure 37, have been obtained by multiplying the predicted 
velocities by a scale factor of 10. The velocities at the computational inlet plane consist, 
in general, of the prescribed incoming wake velocity perturbation and a velocity associated 
with the unsteady pressure response of the FEGV. In this case, although the 1 BPF pressure 
response attenuates with increasing axial distance from the blade row, it still provides a small 
contribution to the unsteady velocity at the inlet plane. The unsteady velocity field, depicted^ 
in Figure 37, shows a strong response at the leading-edge of each blade. Downstream of the 
blade leading edges, the velocity amplitudes diminish with axial distance, and the unsteady 
flow is dominated by counter-rotating, vortical eddies. 

The excitation at 1 BPF produces attenuating acoustic responses upstream and down- 
stream of the FEGV. Thus, the pressure variations near the computational inlet and exit 
boundaries, depicted in Figure 38, are due to acoustic responses that die out with increas- 
ing axial distance from the blade row. The amplitudes of the least-damped 1,0 and 0,0 
acoustic response modes, at the computational inlet and exit boundaries, are 0.00372 and 
0.00242, respectively. The next strongest mode at inlet is the 0,0 mode, with an amplitude 
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of 0.000583. The next strongest mode at exit is the 0,2 mode, which has an amplitude of 
0.000438. However, the exit 0,2 response is more highly damped than the exit 0,1 response. 
The average steady pressure at the FEGV inlet plane is 2.325. Thus, the unsteady pressure 
amplitude of 0.00372 is 0.16% of the average inlet pressure. 

The unsteady blade-surface pressure distributions resulting from the 1 BPF wake exci- 
tation are shown in Figure 39. The surface pressure responses are given for the modified 
FEGV blades, which extend from the blade leading edge near £ = 0 to the sharp trailing 
edge near £ = 1.1 lc^. The trailing edge of the original (blunt-edged) blade occurs near 
£ = c ax . Since the amplitude and phase of the 1 BPF wake excitation vary significantly from 
hub to tip, cf. Figure 32, and this excitation occurs at a moderate reduced frequency, the 
surface-pressure response exhibits a relatively complicated spanwise behavior. A disconcert- 
ing feature of the results, given in Figure 39, are the large pressure gradients predicted in the 
tip-region on the blade suction surface, near the trailing edge. This behavior could be due 
to numerical losses in the steady solution resulting from the prescription of high tip inflow 
angles. 


2 BPF Wake Excitation 

The eigensolutions for the 2 BPF excitation, cf. Figure 40, indicate that four modal 
acoustic responses, the 1,0, 1,1, 1,2 and 1,3 responses, propagate upstream and downstream 
of the FEGV. The axial wave numbers of the 1,0, 1,1, 1,2 and 1,3 upstream propagating 
modes at inlet are 7.884, 7.347, 6.705 and 5.909, respectively. Those of the 1,0, 1,1 1,2 and 
1,3 downstream propagating modes at exit are -2.434, -2.289, -1.736 and -0.645. Of the 
captured modes, the least-damped upstream and downstream attenuating responses have 
axial eigenvalues \ 2 ,o = 3.992 + 3.417* and Xi,a = -1-548 + 2.009*. An attenuating 1,4 mode 
was not determined by the upstream far-field analysis. Because of the higher mean flow 
speeds, there are more propagating modes in response to the 2 BPF excitation at the high 
power setting, than at the low power setting. Also, the axial wave number of the response 
in a given propagating mode has a higher magnitude. 

The imaginary components of the unsteady velocity and pressure resulting from the 2 
BPF excitation are shown in Figures 41 and 42. A scale factor of 25 was applied to the 
predicted velocity, to provide the results given in Figure 41. The velocity distribution along 
the computational inflow boundary consists of the prescribed wake-excitation velocity, along 
with a contribution from the velocity associated with the propagating pressure or acoustic 
response of the FEGV. Again, there are large velocity responses at the leading edge of 
each vane, particularly near the tip, and the unsteady flow becomes dominated by counter- j 
rotating, vortical eddies as it proceeds through the blade row. 

The unsteady pressure responses, depicted in Figure 42, to the 2 BPF wake excitation 
clearly indicates propagating acoustic responses in the far field of the FEGV. The amplitudes 
of the upstream propagating 1,0, 1,1, 1,2 and 1,3 acoustic response modes at inlet are 0.00180, 
0.000869, 0.00231 and 0.000756. The least-damped captured mode, i.e., the 2,0 mode, has an 
amplitude of 0.000302 at the computational inflow boundary. The strongest damped mode at 
inlet is the 0,1 mode, which has an amplitude of 0.000444. At exit, the propagating 1,0, 1,1, 
1,2 and 1,3 modes have amplitudes of 0.00281, 0.00223, 0.00246 and 0.00149. The strongest, 
and least-damped, mode at exit is the 1,4 mode, which has an amplitude of 0.000132. Here, 
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the exit 1,0 mode occurs at the highest amplitude, and this is 0.12% of the average inlet 
pressure. There is no clearly dominant acoustic response mode in either the upstream or 
the downstream regions of the flow. In addition, the modes that contribute to the pressure 
responses at inlet and exit have large radial variations and reach maximum amplitudes at 
different radial locations. Such features make it difficult to analyze the acoustic responses of 
the FEGV, based only on the pressure contours plotted in Figure 42. 

The unsteady surface-pressure responses of the FEGV to the 2 BPF wake excitation are 
shown in Figure 43. The real and imaginary, surface-pressure responses show very com- 
plicated radial behaviors. They also show spurious chordwise variations, along the suction 
surface, near the tip and near the trailing edge, which are indicative of error in the unsteady 
solution. Again, this spurious behavior may be due to the numerical losses associated with 
the high, tip, inflow angles prescribed for the TURBO steady-flow calculation. 

The radial shapes of the propagating acoustic response modes for the 2 BPF excitation 
are shown in Figures 44 and 45. The propagating response disturbances occur in the m = 1 
circumferential mode, i.e., at m = -2 N B + N v = 9 or o x = 2nfh/N v - 72 deg. Hence, they 
repeat nine times around the wheel and travel circumferentially counter to the direction 
of rotor rotation. Each of the propagating modes, including those at // = 0, show strong 
spanwise pressure variations. 

The predicted sound pressure levels, in dB, for the 1,0, 1,1, 1,2 and 1,3 acoustic responses 
at inlet are 124.8, 119.2, 127.2 and 117.4, respectively. Those for the 1,0, 1,1, 1,2 and 1,3 
acoustic responses at exit are 127.4, 126.5, 126.6 and 122.4. The total response sound pressure 
levels, at inlet and exit are 129.9 dB and 132.1 dB, respectively. The predicted modal and 
total sound power levels of the acoustic responses at inlet and exit, and, in parentheses, the 
measured [Hei99] sound power levels of the responses at exit, are given in Figures 44 and 45. 
The agreement between the analytical and the experimental results for the exit sound power 
levels is excellent. However, considering the various approximations made at this stage of the 
development of the LINFLUX Prediction System, this level of agreement must be regarded 
as somewhat fortuitous. Again, acoustic responses at small amplitudes produce significant 
sound power and sound pressure levels. At the high power setting, 2 BPF acoustic response 
amplitudes at approximately 0.1% of the average inlet pressure, produce sound pressure and 
sound power levels on the order of 120 dB. 

Animations 

We have constructed animations to illustrate the time-dependent behaviors of the LIN- 
FLUX solutions for the 1 and 2 BPF rotor-wake/FEGV interactions, at the high ADP power 
setting. The animations include a 3D view of the unsteady blade-surface pressures over the-3 
entire FEGV, and 2D radial views of the axial velocity and pressure at the 94% spanwise 
station, over 4 blade passages. 

For the 1 BPF excitation, the 3D animation shows an unsteady surface-pressure response, 
that travels from blade to blade, counter to the direction of rotor rotation. Pressure distur- 
bances at high amplitude occur on the suction surfaces of the blades, near the tip, and travel 
upstream along the blade surface. On the pressure surface, a high-amplitude response occurs 
in the middle of the blade, travels outward to the tip, and then upstream to the leading edge. 

The 2D axial-velocity animation shows well-defined rotor-wakes that extend from the 
computational inlet plane to the FEGV leading-edge plane. The wake excitation repeats 18 
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times around the wheel, and travels circumferentially in the direction of rotor rotation and 
axially towards the FEGV. Downstream of the FEGV leading-edge plane, the axial velocity 
perturbations have a very complicated behavior, which depends on the interaction between 
the prescribed excitation and the FEGV. Within the blade passages, these disturbances show 
strong blade-to-blade motions that are associated with the unsteady pressure response of 
the FEGV. Aft of the leading-edge plane, the axial velocity disturbances travel downstream, 
essentially, in the axial flow direction. 

The 1 BPF pressure animation indicates high amplitude pressure variations within each 
blade passage. These emanate from the blade suction surfaces, and travel into the passage 
and upstream. Upstream of the FEGV, the pressure disturbances travel counter to the direc- 
tion of rotor rotation and away from the blade row. Downstream, the pressure disturbances 
travel in the direction of rotor rotation and, seemingly, towards the blade row. In both 
regions, the pressure disturbances attenuate with increasing axial distance from the blading. 

For the 2 BPF excitation, the 3D animation shows an unsteady surface-pressure response 
that travels circumferentially, counter to the direction of rotor rotation. On both the blade 
suction and pressure surfaces, the pressure disturbances travel from trailing to leading edge 
and from hub to tip. 

The 2D, axial velocity animation indicates well-defined rotor wakes from the inlet plane 
to the vane leading-edge plane. The wake excitation repeats 36 times around the wheel, and 
travels in the direction of rotor rotation and towards the blade row. This excitation breaks 
down aft of the vane leading-edge plane, where the axial velocity perturbations become 
weaker and less organized, and travel downstream, essentially in the axial direction. The 2 
BPF, time-dependent, axial-velocity perturbation indicates a much smaller blade-to-blade 
movement, than that indicated for the 1 BPF excitation. 

The 2 BPF, unsteady pressure animation shows strong pressure activity at the blade 
surfaces and a complicated pressure behavior within the blade passages, where the pressure 
disturbances generally travel upstream. The behavior of the pressure perturbations changes 
rather abruptly to that typical of far-field pressure responses just upstream of the FEGV 
leading-edge plane and just downstream of the trailing-edge plane. The outward propagating 
pressure responses upstream and downstream of the blade row are strong and well defined. 
They travel circumferentially counter to the direction of rotor rotation, and axially, away 
from the blade row. The direction of propagation of the upstream disturbance is at a 
small angle, measured clockwise, from the negative axial direction, that of the downstream 
disturbance, is at an angle of roughly 45 deg, measured counter-clockwise, from the positive 
axial direction. 
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9. Concluding Remarks 


The LINFLUX analysis is being developed to serve as the unsteady aerodynamic compo- 
nent of advanced, aeroacoustic and aeroelastic, design prediction systems for turbomachines. 
It is based on the linearized Euler equations, an implicit, wave-split, finite- volume analysis to 
determine the unsteady flow in the near field of a blade row, and eigenanalyses to determine 
the unsteady flows in the far field. The near- and far- field solutions are coupled at the inflow 
and outflow boundaries of the near-field domain. To date, the theoretical effort has focused 
on developing the unsteady aerodynamic model and the near- and far-field solution proce- 
dures; implementing these into two- and three-dimensional, unsteady-flow codes; validating 
these codes for benchmark unsteady flows; and, under the present effort, demonstrating the 
3D code, for aeroacoustic design, via application to realistic wake/blade-row interactions. 

As noted previously, the application of LINFLUX to predict the unsteady aerodynamic 
response of a blade row to a wake excitation involves the use of four codes: TIGER, to 
generate a 3D H-grid; AWAKEN, to define steady and unsteady inflow conditions; TURBO, 
to determine the steady background flow; and LINFLUX to determine the unsteady flow. 
Under the present effort, we have applied this system of codes to predict steady and unsteady 
flows through the exit guide vane of a NASA/PW advanced ducted propulsor. The unsteady 
flows are excited by circumferential nonuniformities in the flow exiting the fan rotor, which 
are due primarily to the rotor wakes. We have considered two ADP power settings, a low 
power or approach setting and a high power or take-off setting, and have applied LINFLUX 
to predict the aerodynamic responses of the FEGV to Fourier components, at multiples of 
the blade passing frequency, of the rotor wake excitations. The latter are based on velocity 
measurements, taken, in an axial plane, downstream of the rotor, at NASA Glenn Research 
Center. 

A number of liberties were taken in modeling the actual rotor- wake/FEGV interac- 
tion, including geometric changes to the vanes and the aft duct, and modifications to the 
circumferentially-averaged, measured, rotor-exit flow. The major approximations are those 
associated with modeling the wake excitation. The circumferential wake velocity distribu- 
tions have been modified in an attempt to remove endwall phenomena and provide a more 
or less classical wake excitation from hub to tip. In addition, empirical correlations have 
been applied to estimate the strength and circumferential distribution of the viscous wake 
excitation at the FEGV midspan leading-edge plane, based on the velocity measurements 
taken upstream. The estimated velocity distribution at this reference axial plane is used to 
determine the complex amplitudes of the Fourier wake excitations that are input, at inlet, _ 
to the inviscid LINFLUX code. Finally, a quasi-3D or strip-theory approximation has been 
employed, in which, at each radial station, the inlet wake velocity perturbation is assumed to 
be an axially and circumferentially convected disturbance that carries no pressure or density. 

We have presented numerical results that describe the flows through the FEGV. These 
include AWAKEN results for the steady inflow conditions and the wake-excitation velocity; 
TURBO predictions for the steady background flow; and LINFLUX predictions for the near- 
field unsteady velocity and pressure, the modal acoustic responses at inlet and exit, and the 
modal and total, sound pressure and sound power, levels in the far field. Unsteady response 
predictions, for the low power setting, have been determined for rotor-wake excitations at 
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1, 2, and 3 BPF. However, the 3 BPF results are not accurate. Unsteady predictions, for 
the high power setting, have been determined for wake excitations at 1 and 2 BPF. In this 
case, a converged LINFLUX solution could not be determined for the 3 BPF excitation. The 
difficulties at 3 BPF are believed to be due to the errors associated with representing the 
highly radially dependent 3 BPF wake disturbances in a quasi-3D manner. 

The steady background flows and the Fourier wake excitations vary significantly with 
radius. Also, the 1 and 2 BPF excitations occur at sufficiently high frequencies that they 
produce propagating or lightly attenuated acoustic responses that are highly dependent on 
radius. These features make it difficult to assess the LINFLUX solutions, particularly in the 
absence of any previous analytical results or detailed near-field response data. The results for 
the 1 and 2 BPF excitations, at both power settings, indicate complicated three-dimensional, 
fluid-dynamic, response phenomena. However, for the most part, these results appear to be 
very reasonable. 

The LINFLUX predictions, for the two power settings, indicate that the 1 BPF wake 
excitations produce attenuating acoustic responses, upstream and downstream of the FEGV 
and, hence, no sound far from the blade row. The 2 BPF excitations produce propagating 
acoustic responses, in both far-field regions, that carry energy away from the blade row. The 
2 BPF excitations occur at similar reduced frequencies. However, the steady flow speeds are 
higher and the wake excitations are stronger at the high power setting. As a result, more 
propagating response modes occur at this setting. These travel away from the blade row at 
higher axial wave numbers, and they produce higher sound levels. At the low power setting, 
three modal responses propagate upstream and two propagate downstream. At the high 
setting, four modal responses propagate upstream and four propagate downstream. The 
propagating acoustic responses at 2 BPF repeat nine times around the wheel, and travel 
counter to the direction of the rotor rotation and away from the blade row. 

The propagating modes all have strong radial variations and significant sound power 
levels. In particular, at the low power setting, the 2 BPF excitation produces an acoustic 
response at exit having a predicted total sound power level of 111.8 dB. The sound power 
levels of the 1,0 and 1,1 modal responses are 110.5 and 106.0 dB. The corresponding measured 
values are 114.8, 110.4 and 112.8 dB, respectively. At the high power setting, the 2 BPF 
excitation produces a predicted total sound power level of 125 dB at exit. The predicted 
sound power levels for the exit 1,0, 1,1, 1,2 and 1,3 propagating acoustic responses are 
121.1, 120.4, 119.2 and 112.5 dB respectively. The corresponding measured power levels are 
126.0, 122.1, 121.4, 119.4 and 110.1 dB, respectively. The 2 BPF excitations produce acoustic 
response amplitudes that are on the order of only 0.1% of the average inlet pressure, but such 
responses produce acoustic power levels on the order of 110 to 120 dB. With the exception 
of those for the 1,1 mode at the low power setting, the predicted and measured modal 
sound power levels at exit are in excellent agreement. However, although very encouraging, 
this level of agreement must be regarded as somewhat fortuitious, because of the various 
approximations made in applying LINFLUX, particularly those associated with modeling 
the wake excitation. 

Animations of the the LINFLUX axial-velocity solutions, for the the 1 and 2 BPF wake 
excitations at the high power setting, indicate well defined wake excitations extending from 
the computational inflow boundary to just upstream of the FEGV midspan leading-edge 
plane. Downstream of this axial plane, the wake excitations break down and the axial 
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velocity perturbations exhibit very complicated behaviors. Within the blade passages, the 
axial- velocity disturbances show strong blade-to-blade motions. Aft of the blade row, weaker, 
less-organized, axial velocity perturbations travel downstream, essentially in the mean-flow 
(axial) direction. The 1 BPF pressure animation shows attenuating pressure responses that 
travel counter to the direction of rotor rotation, upstream of the FEGV, and in the direction 
of rotor rotation, downstream. The 2 BPF pressure animation indicates strong, well-defined, 
propagating pressure responses upstream and downstream of the FEGV. These travel counter 
to the direction of rotor rotation and away from the blade row. 

Based upon the validation studies conducted to date and the present application to fan- 
wake/FEGV interactions, it appears that the 3D LINFLUX analysis will be a useful predic- 
tion capability for determining the unsteady pressure responses of blade rows to the various 
sources of unsteady excitation. Therefore, it should be a valuable resource for aeroacoustic 
design, and for understanding the flow physics associated with blade-row noise generation. 
The results given in this report are based on an early application of the LINFLUX Predic- 
tion System to realistic wake/blade-row interactions. It is anticipated that solutions will 
improve as experience is gained in using the various codes. More importantly, with further 
development, the computational resources required by the TURBO and LINFLUX analyses 
will decrease, and various aspects of the unsteady flow model will be improved. 

The major limitations, associated with applying the present version of LINFLUX to 
wake/blade-row interactions, are those related to the wake excitation model. At present, 
hub and duct phenomena are ignored, the effects of viscous diffusion are modeled empir- 
ically, and the wake excitation, at inlet, is represented as a quasi-3D, convected, velocity 
disturbance. Future work should be focused on improving these aspects of the wake/blade- 
row interaction analysis. Improved representations of the excitation could be achieved by 
having sufficient experimental or computational data on the rotor-exit flow to accurately 
specify the thermodynamic properties and the axial behavior of the wake excitation at the 
stator inlet. In addition, although it would entail a major effort, a viscous version of LIN- 
FLUX would enable the effects of wake diffusion to be accommodated within the unsteady 
aerodynamic analysis. 

Several computational issues should also be investigated to improve the performance of 
the LINFLUX code. These include improvements to the far-field eigenanalysis to ensure that 
all important acoustic response modes are captured. It would also be of value to develop a 
far-field analysis for variable radii ducts. Finally, the development of a parallel version of 
LINFLUX would allow the grids, needed to resolve unsteady flows at high subsonic Mach 
numbers and/or high excitation frequencies, to be readily accommodated. 
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Figure 1: Rotating axial compressor blade row operating within an annular duct. 
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Figure 2: Nomenclature for wake/blade-row interaction. 
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Figure 6: FEGV computational grid at 4 radial stations. 
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Figure 7: Steady flow properties far upstream of the FEGV operating at the low power 
(approach) setting. 
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Figure 8: Absolute and relative steady velocities at the rotor exit (FEGV inlet) for the low 
power setting. - 


NASA/CR— 2001-210713 


52 















30 . 



- 10 . 



Figure 10: Fourier components of the analytic, rotor-wake excitation velocity, 5 R , for the 
low power setting. 
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Figure 12: Steady surface-pressure distributions at four radial stations for the FEGV oper 
ating at the low power setting. 
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Figure 13: Steady flow properties at the computational inlet and exit planes for the FEGV 
operating at the low power setting. 
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Figure 16: Imaginary component of the unsteady pressure at four radial stations for the 
FEGV operating at the low power setting and subjected to the 1 BPF wake excitation. 
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Figure 17: Unsteady surface-pressure distributions at four radial stations for the FEGV 
operating at the low power setting and subjected to the 1 BPF wake excitation. 
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Figure 18: Axial eigenvalues, x = P + of acoustic disturbances far upstream and far 
downstream of the FEGV operating at the low power (approach) setting and subjected to 
an unsteady excitation at 2 BPF (u = 7.406 and a = —288.0 deg). 
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Figure 19: Imaginary component of unsteady velocity field at four radial stations (scale 
factor = 30) for the FEGV operating at the low power setting and subjected to the 2 BPF 
wake excitation. " ‘ 
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Figure 20: Imaginary component of unsteady pressure field at four radial stations for the 
FEGV operating at the low power setting and subjected to the 2 BPF wake excitation. 
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Figure 21: Unsteady surface-pressure distributions at four radial stations for the FEGV 
operating at the low power setting and subjected to the 2 BPF wake excitation. 


NASA/CR— 2001-210713 


64 





Inlet: PWL tot = 106.1 db 



v / r o r/rxj r/ro 


Exit: PWL tot = 111.8 (114.8) db 



Figure 22: Radial pressure modes and modal and total sound power levels of the propagating 
acoustic responses in the far field of the FEGV operating at the low power (approach) setting 
and subjected to the 2 BPF (u; = 7.406 and a = —288.0 deg) wake excitation. 
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Figure 23: Axial eigenvalues, x = P + of acoustic disturbances far upstream and far 
downstream of the FEGV operating at the low power (approach) setting and subjected to 
an unsteady excitation at 3 BPF ( u — 11.109 and a = —432.0 deg). 
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Figure 24: Imaginary component of the unsteady velocity at four radial stations (scale factor 
= 75) for the FEGV operating at the low power setting and subjected to the 3 BPF wake 
excitation. 
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Figure 25: Imaginary component of the unsteady pressure at four radial stations for the 
FEGV operating at the low power setting and subjected to the 3 BPF wake excitation. 


NASA/CR— 2001-210713 


68 



0.0030 


0.0015 
Re {p} 

0.0000 

-0.0015 

-0.0030 
0.0030 

0.0015 
Im {p} 

0.0000 

-0.0015 

-0.0030 

Figure 26: Unsteady surface-pressure distributions at four radial stations for the FEGV 
operating at the low power setting and subjected to the 3 BPF wake excitation. 
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Figure 27: Radial pressure modes and modal and total sound power levels for the propagating 
acoustic responses far upstream of the FEGV operating at the low power (approach) setting* 
and subjected to the 3 BPF (u = 11.109 and a - -432.0 deg) wake excitation. 
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Figure 28: Radial pressure modes and modal and total sound power levels for the propagating 
acoustic responses far downstream of FEGV operating at the low power (approach) setting 
and subjected to the 3 BPF (a; = 11.109 and o = — 432.0deg) wake excitation. 
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Figure 29: Steady flow properties far upstream of the FEGV operating at the high power 
(take-off) setting. 
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Figure 30: Absolute and relative steady velocity 
power setting. 
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Figure 31: Analytic, rotor- wake, tangential velocity perturbation at the FEGV midspan 1 
leading-edge plane (£ = 0) for the high power (take-off) setting. 
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Figure 32: Fourier components of the analytic, rotor-wake excitation velocity, v R , for the 
high power setting. 
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Figure 34: Steady surface-pressure distributions at four radial stations for the FEGV oper 
ating at the high power setting. 
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Figure 35: Steady flow properties at the computational inlet and exit planes of the FEGV 
operating at the high power setting. 
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Figure 36: Axial eigenvalues, \ — P + > of acoustic disturbances far upstream and far 

downstream of the FEGV operating at the high power (take-off) setting and subjected to 
an unsteady excitation at 1 BPF ( uj = 3.677 and a = —144.0 deg). 
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Figure 37: Imaginary component of 
= 10) for the FEGV operating at t' 
excitation. 
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Figure 38: Imaginary component of the unsteady pressure at four radial stations for the 
FEGV operating at the high power setting and subjected to the 1 BPF wake excitation. 
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Figure 39: Unsteady surface-pressure distributions at four radial stations for the FEGV 
operating at the high power setting and subjected to the 1 BPF wake excitation. 
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Figure 40: Axial eigenvalues, x = /? + of acoustic disturbances far upstream and far 
downstream of the FEGV operating at the high power (take-off) setting and subjected to 
an unsteady excitation at 2 BPF (u = 7.355 and o = —288.0 deg). 
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Figure 41: Imaginary component of the unsteady velocity at four radial stations (scale factor 
= 25) for the FEGV operating at the high power setting and subjected to the 2 BPF wake 
excitation. 
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Figure 43: Unsteady surface- pressure distributions at four radial stations for the FEGV 
operating at the high power setting and subjected to the 2 BPF wake excitation. 
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Figure 44: Radial pressure modes and modal and total sound power levels for the propagating 
acoustic responses far upstream of the FEGV operating at the high power (take-off) setting 
and subjected to the 2 BPF (w = 7.355 and a = —288.0 deg) wake excitation. 
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Figure 45: Radial pressure modes and modal and total sound power levels for the propagating 
acoustic responses far downstream of the FEGV operating at the high power (take-off) setting 
and subjected to the 2 BPF (u = 7.355 and a = —288.0 deg) wake excitation. 
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